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THEORETICAL  AND  EXPERIMENTAL  STUDY  OF  THERMOACOUSTIC  ENGINES 


1.0  Introduction 

The  major  thrust  of  this  research  is  the  use  of  ceramic  elements  as  the  stack 
structure  in  a  thermoacoustic  device  and  the  development  of  an  acoustics  based  theory  to 
calculate  the  behavior  of  the  thermoacoustic  engine.  This  work  builds  on  previous  work 
performed  at  the  Naval  Postgraduate  School  and  Los  Alamos  National  Laboratory  on 
thermoacoustic  engines  and  refrigerators  as  well  as  the  Physical  Acoustics  Research  Group 
extensive  work  on  sound  propagation  through  porous  media. 

The  main  body  of  this  report  describes  the  general  formulation  of  thermoacoustics 
for  stacks  having  arbitrary-shaped  pore  cross-sections.  Attached  appendices  are  papers 
describing  the  work  performed  by  Dr.  Bass  at  the  Naval  Postgraduate  School,  the 
development  and  testing  of  the  theory  of  sound  propagation  in  square  pore  material,  and  the 
design  parameter  of  the  UM  low  frequency  thermoacoustic  driver. 
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2.0  GENERAL  FORMULATION  OF  THERMOACOUSTICS  FOR  STACKS 
HAVING  ARBITRARILY-SHAPED  PORE  CROSS-SECTIONS 

2.1  Introduction 

The  subject  of  modem  thermoacoustics  is  primarily  due  to  the  work  of  N.  RotL  His 
review  paper^  describes  much  of  this  work.  Thermoacoustic  streaming  in  a  driven 
resonance  tube  was  also  investigated  by  Merkli  and  Thomann.^  Rott  and  Merkli  and 
Thomann  were  interested  in  thermoacoustic  effects  in  a  single  tube  having  a  circular  cross 
section.  Rott  and  Zouzoulas^  also  investigated  thermally  driven  acoustic  oscillations  for 
circular  tubes  with  variable  cross-sectional  area. 

Wheatley,  Cox,  Swift,  Hoffler  and  others  have  approached  thermoacoustics  from  a 
more  traditional  thermodynamics  point  of  view.  Much  of  this  work  is  reviewed  by  Swift.'^ 
They  use  thermoacoustic  elements  such  as  heat  exchangers  and  a  stack  as  shown  in  Fig. 
2.1a  to  investigate  thermoacoustically-driven  prime  movers  and  refrigerators.  An  exposed 
view  of  a  thermoacoustic  element  is  shown  in  Fig.  2.1b.  Thermoacoustic  elements  consist 
of  a  parallel  combination  of  many  elementary  capillary  tubes  or  pores.  For  example  in  Fig. 
2.1b  the  pores  have  square  cross-sections.  The  theory  for  a  thermoacoustic  heat  engine  is 
built  up  from  knowledge  of  thermoacoustic  effects  in  a  single  capillary  tube.  Both  Rott^ 
and  Swifti^  considered  thermoacoustic  effects  for  acoustic  oscillation  between  parallel 
plates.  Merkli  and  Thomann^  and  Rott^-^  investigated  thermoacoustic  effects  for  tubes 
with  circular  cross-sections. 

Our  initial  intent  was  to  investigate  thermoacoustic  effects  in  a  stack  having  square 
pores  as  shown  in  Fig.  2.1b.  An  inexpensive  source  of  square  pore  stack  material  is  the 
ceramic  plug  used  in  automobile  catalytic  converters.^  We  previously  used  this  ideal- 
geometry  material  to  investigate  sound  propagation  in  porous  media.^  Since  ceramic  has  a 
low  thermal  conductivity  in  comparison  to  metals,  it  is  ideal  for  use  as  a  stack  where  heat 
loss  due  to  normal  conduction  down  the  stack  reduces  efficiency. 

In  this  paper,  we  investigate  thermoacoustics  for  stacks  having  arbitrary  pore 
geometries,  (parallel  slits,  rectangular  pores,  equilateral-triangle  pores,  circular  pores,  etc). 
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In  particular,  we  are  interested  in  the  following  question:  What  are  the  minimum  necessary 
calculations  to  describe  the  acoustics  of  gas-filled  thermoacoustic  elements  made  of 
arbitrary-perimeter  capillary  tubes.  An  example  of  an  arbitrary-perimeter  capillary  tube  is 
shown  in  Fig.  2.1c.  Rott^  pursued  this  question  to  the  point  of  computing  the  acoustic 
field  quantities  pressure,  density,  particle  velocity,  and  temperature  for  arbitrary  perimeter 
tubes,  and  gave  examples  of  parallel  slit  and  circular  pore  geometries.  Here  we  consider 
the  acoustic  field  quantities  and  the  second  order  energy  flow  for  arbitrary  perimeter  pores. 
Heat  and  work  flow  are  compared  in  the  short  stack  approximation  for  stacks  having  the 
aforementioned  pore  geometries.  In  addition,  we  make  connections  between  thermo¬ 
acoustic  theory  and  capillary-tube-based  porous  media  theory. 

Once  the  acoustical  behavior  of  the  various  thermoacoustic  elements  has  been  worked 
out,  the  elements  must  be  connected  in  series  inside  of  a  resonator  as  shown  in  Fig.  2. 1  a. 
Previously,  numerical  integration  of  the  acoustical  equations  was  used  to  compute  field 
quantities  in  the  stack  since  in  general  a  temperature  gradient  exists  from  one  side  to  the 
other.^  The  physical  parameters  of  ambient  density,  viscosity,  specific  heat,  sound  speed, 
thermal  conductivity,  etc,  are  temperature  dependent  and  thus  depend  on  location  within  the 
stack.  An  alternative  method  developed  here  assumes  that  the  stack  may  be  broken  up  into 
isothermal  subsections  over  which  the  physical  parameters  are  constant.  When  this 
assumption  is  made,  techniques  available  to  analyze  waves  in  layered  media  are  applicable. 
In  particular,  specific  acoustic  impedance  and  pressure  translation  theorems  are  developed 
to  compute  all  acoustical  field  quantities  and  energy  flow  at  each  point  in  the  resonator 
shown  in  Fig.  2.1a.  This  method  of  solution  is  easily  generalized  to  more  complicated 
systems,  e.g.  a  resonator  containing  a  refrigeration  stack,  a  prime  mover  stack,  and  4  heat 
exchangers.  These  translation  theorems  are  uniformly  applicable  to  propagation  in  heat 
exchangers,  stacks,  and  in  the  wide-tube  sections  of  the  resonator. 
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acoustic  heat  heat 

driver  exch  stack  exch 


resonance  tube 

a) 


exposed  view  of  a  thermoacoustic  element 


Tambient  =  To  +  TozZ 


single  arbitrary-perimeter  tube  of  a 
thermoacoustic  element 


Figure  2.1  (a)  Generic  arrangement  used  in  thermoacoustics.  Thermoacoustics 
elements  are  the  heat  exchanges  and  stack,  (b)  An  exposed  view  of 
a  thermoacoustic  element  consisting  of  a  parallel  combination  of 
square  capillary  tubes,  (c)  A  single  arbitrary-perimeter  capillary 
tube  for  use  in  a  thermoacoustic  element. 
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2.2  Propagation  in  Thermoacoustic  Elements  from  a  Porous  Media  Analysis 

A  common  design  for  the  stack  and  heat  exchangers  in  a  thermoacoustic  engine  are 
equally  spaced  parallel  plates.  Fluid  undergoes  oscillitory  flow  between  the  plates  allowing 
for  greater  heat  transfer  than  would  occur  if  the  fluid  was  static.  The  pores  in  the  parallel 
plate  arrangement  are  described  in  their  transverse  direction  as  parallel  slits  and  as  straight 
tubes  in  their  longitudinal  direction.  For  arbitrary  geometries,  the  stack  or  heat  exchangers 
can  be  described  as  a  section  of  a  porous  media.  In  this  section,  we  first  establish  the  fluid 
field  equations  and  assumptions  necessary  to  treat  the  general  case.  A  wave  equation  for 
the  pressure  in  a  single  pore  is  established.  Enroute,  we  make  reference  to  the  terminology 
used  in  acoustical  modeling  of  porous  media.  With  the  specific  acoustic  impedance 
assumed  known-  at  the  hot  end  of  the  stack,  we  derive  impedance  and  pressure  translation 
theorems  for  the  stack.  Heat  and  work  flow  are  computed  for  arbitrary  pore  geometries 
and  are  expressed  in  terms  of  pressure  and  specific  acoustic  impedance. 

A.  Fluid  field  equations  and  assumptions. 

The  transverse  coordinates  in  a  pore  are  taken  to  be  x  and  y,  and  the  longitudinal 
coordinate  is  z.  For  the  stack,  a  positive  linear  temperature  gradient  exists  on  the  pore  in 
the  z  direction.  Heat  flows  to  and  from  the  gas  inside  the  pore  to  the  walls  so  that,  on  the 
average,  the  temperature  along  the  pore  is  the  same  for  the  fluid  and  the  pore  walls  at  the 
boundary.  We  assume  that  the  pore  walls  are  of  sufficiently  high  heat  capacity  that  the  pore 
wall  temperature  is  unaffected  by  temperature  variations  in  the  fluid  caused  by  an  acoustic 
wave.  We  also  assume  constant  frequency  pressure  variations  exist  in  the  pore  and  that  the 
pore  walls  are  rigid.  The  pore  is  taken  to  be  infinitely  long  in  the  z  direction.  With  these 
assumptions,  we  would  like  to  completely  specify  the  acoustic  field  in  the  pore  to  first 
order  in  the  acoustic  variables. 


5 


The  fluid  quantities  in  a  pore  approximated  to  first  order  are: 


P(z,t)  =  Po  +  Pi(z)  exp{-\m) ,  (2.1) 

v(x,y,z,t)  =  v(x,y,z)  ej:p(-icot)  ,  (2.2) 

T(x,y,z,t)  =  To  ^  exp(-\m) 

=  To  +  Toz  z  +  Ti(x,y,z)  ej:p(-icot),  (2.3) 

s(x,y,z,t)  =  so(z)  +  si(x.y,z)  exp{-\m) ,  (2.4) 

p(x,y,z,t)  =  po(z)  +  pi(x,y,z)  exp(-icDt)  .  (2.5) 


Because  of  the  ambient  temperature  gradient,  physical  parameters  including  density, 
viscosity,  thermal  conductivity,  adiabatic  sound  speed,  coefficient  of  thermal  expansion, 
and  the  specific  heats  also  depend  on  position.  We  ignore  these  variations  within  a  section 
of  small  length.  Later  we  will  include  this  temperature  dependence  by  treating  the  stack  as 
a  sum  of  short  stacks  with  different  temperatures.  For  physical  parameters  independent  of 
position,  the  fluid  field  equations  to  first  order  are^ 

-i  CO  po  vz(x,y,z)  =  -  +  Tj  vz(x,y,z)  .  (2.6) 


i  CO  pi(x,y,z)  ^ ^  p,(,)  ^  ,  o  ^2.7) 


pl(x,y,z)  =  -po(z)  pTi(x,y,z)  +  XPi(z)  , 

c^ 


(2.8) 


si(x,y,z) 


^Ti(x,y,z)  -  —  Pi(z)  , 
PO 


(2.9) 
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and 


-  i(Opo(z)CpTi(x,y,z)  +  po(z)CpVz(x,y,z)Toz  =  -i(Dl3ToPi(z)  +  KV2xyTi(x,y,z),  (2. 10) 

where  the  transverse  Laplacian  operator  is  defined  by  =  (82/3x2  +  82/8y2).  in 
order,  these  equations  approximately  express  the  z  component  of  the  equation  of  motion, 
continuity  or  mass  conservation,  equations  of  state  for  density  and  entropy,  and  heat 
transfer.  Except  for  the  Tqz  terms  in  Eq.  (2.10),  these  are  the  equations  for  the  low 
reduced  frequency  solution^  given  by  Zwikker  and  Kosten^O  in  their  solution  for  the 
propagation  of  sound  in  circular  pores.  In  the  derivation  of  Eq.  (2.10),  the  convective 
derivative  for  the  entropy  is  evaluated  using  the  equation  of  state  (2.9)  and  the  relation  v- 
Vs  =  Vz(x,y,z)  dso(z)/dz  =  Cp  Vz(x,y,z)  Tqz  /Tq.  Equation  (2.10)  shows  that  the 
temperature  at  a  fixed  position  changes  due  to  motion  of  the  ambient  fluid,  due  to 
compression  of  the  gas,  and  due  to  heat  conduction. 

In  writing  Eq.  (2.6),  we  assume  that  the  longitudinal  fluid  velocity  is  much  larger 
than  the  transverse  velocity.  We  also  assume  that  the  transverse  variation  of  the 
longitudinal  velocity  is  much  greater  than  the  longitudinal  variation.  The  latter  assumption 
is  sensible  since  the  acoustic  wavelength  is  much  larger  than  the  pore  diameter.  The  fluid 
velocity  vanishes  at  the  pore  wall  and  changes  rapidly  in  a  viscous  boundary  layer  away 
from  the  pore  wall.  In  Eq.  (2.10),  we  have  also  made  the  assumption  that  the  transverse 
variations  of  temperature  are  much  greater  than  the  longitudinal  variation.  Further 
discussion  of  these  approximations  may  be  found  in  references  4  and  9. 

Several  variations  in  notation  should  be  noted.  First,  the  viscous  and  thermal 
boundary  layer  thicknesses  are  given  by  5v  =  (2T|/o)po)*/2  and  5k  =  (2K/(Opocp)^/2  Swift^ 
writes  most  equations  in  terms  of  5vand  5k.  Tijdeman^  and  Attenborough^^  introduce  a 
dimensionless  "shear  wave  number"  X  =  R(poo)/T|)fr2  or  A.  =  2^2  R/5y  where  R  is  a 
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characteristic  transverse  dimension  of  the  pore.  They  also  use  the  dimensionless  thermal 
disturbance  number  Xj  =  R(poCi)Cp/K)l/2  qj-  for  the  ratio  of  the  pore  radius 

to  the  thermal  boundary  layer  thickness.  Use  of  the  Prandtl  number  Np,-  =  rjCp/K  gives  the 
relation  A,x  =  >.  Npr^^.  For  definiteness,  we  take  R  to  be  twice  the  ratio  of  the  transverse 
pore  area  to  the  pore  perimeter  so  for  a  circular  or  square  pore,  R  is  just  the  pore  radius.  In 
this  paper  we  use  the  X,  and  X,t  notation. 

B.  Transverse  velocity  profile  in  a  pore. 

To  obtain  a  solution  to  the  equation  of  motion,  Eq.  (2.6),  the  z  component  of  velocity 
is  taken  to  be 


vz(x,y,z) 


F(x,y;X)  dPi(z) 

icopo  ^ 


(2.11) 


From  the  equation  of  motion,  Eq.  (2.6),  F(x,y,X.)  satisfies 


F(x,y;X)  +  ^  V2xy  F(x,y;X)  =  1 


(2.12) 


subject  to  the  boundary  condition  that  F(x,y,A,)  is  zero  ai  the  pore  walls.  As  we  shall  see, 
this  is  the  only  differential  equation  need  be  solved. 

Denote  by  V7(z)  and  F(X.)  the  average  of  Vz(x,y,z)  and  F(x,y;X,)  over  the  cross  section 
of  the  pore.  The  equation  of  motion  for  the  fluid  can  now  be  expressed  simply  as 

vz(z)  =  F(k).  (2.13) 

icopo 
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The  traditional  approach  of  porous  media  modeling^O.ll  introduces  a  complex  density  at 
this  point  which  is  defined  as 


p(z;>.) 


Po(z) 

F(X) 


(2.14) 


Use  of  the  terminology  complex  here  implies  that  p(z;X)  is  a  complex  number.  The 
complex  density  is  the  apparent  dynamical  density  of  the  fluid  in  the  pore. 


C.  Transverse  temperature  profile  in  a  pore:  complex  compressibility. 

The  excess  temperature  in  the  pore  fluid  is  given  by  Eq.  (2.10).  Algebraic 
rearrangement  followed  by  use  of  the  thermodynamic  relation  Tq  /  Cp  =  (y-  1)  /  c^  for 
the  first  term  on  the  right  and  by  use  of  Eq.  (2.11)  results  in  a  simplified  expression  for  the 
excess  temperature, 


Ti(x,y,z)  +  ^V2,y  Ti(x,y.z)  -^Pi(z) -^F(x,y;>.)  .  (2.15) 


c^pop 


pOO)^ 


We  may  assume  a  superposition  form  of  Ti(x,y,z), 
Ti(x,y,z)  =  Ta(x,y)  Pi(z)  +  'tb(x,y)  . 


(2.16) 


Following  Rott*,  Eq.  (2.16)  is  split  into  two  equations  corresponding  to  the 
inhomogeneous  terms, 


'Ca(x,y)  + 


R2 

i 


V^xy  Xa(x,y) 


y-  1 

c2  po  p 


(2.17) 
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(2.18) 


'Cb(x,y) 


R2 
i  A.x^ 


V^xy  Tb(x.y) 


-^F(x,y;>.)  . 


Comparing  Eq.  (2.17)  with  Eq.  (2.8),  the  solution  for  Eq.  (2.17)  follows  immediately. 


ta(x,y)= 


Y-  1 
c2  po  P 


F(x,y;>.T) 


(2.19) 


To  obtain  the  solution  for  Xb(x,y)  in  Eq.  (2.18),  Eq.  (2.19)  can  be  used  in  Eq.  (2.17)  along 
with  the  differential  equation  for  F(x,y;X)  in  Eq.  (2.12)  to  prove  the  following  result. 


'Cb(x,y)  =  - 


Tqz  F(x,y;XT)  -  Npf  F(x,y;>.) 


PO 


1  -  N 


pr 


(2.20) 


The  general  solution  for  the  excess  temperature  from  Eq.  (2.16)  is 


Ti(x,y,z)  =  ^ ^  F(x,y;^T)  Pl(z) 
c2  Po  p 

_  Tqz  F(x,y;>.T)  -  Npr  F(x,y;X)  dPi(z) 

poco2  1  -  Npr  ^ 


The  equation  of  state  for  the  acoustic  density  fluctuation,  Eq.  (2.8),  can  be  combined  with 
the  expression  for  the  acoustic  temperature  ,  Eq.  (2.21)  to  give 


pl(x,y,z)  =  ^  (y  -  (y-  1)  F(x,y,XT))  Pl(z) 

PTqz  F(x,y;XT)  -  Npr  F(x,y;X)  dPi(z) 

(o2  1  -  Npr  ^ 


(2.22) 
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The  first  term  on  the  right  of  Eq.  (2.22)  can  be  used  to  define  a  complex  compressibility. 

When  Toz  is  zero,  Eq.  (2.22)  can  be  used  to  define  a  complex  compressibility  which 
is  useful  in  porous  media  theory.^Odl  Denote  by  pi(z),  F(Xx).  F(A.)  as  the  average  of 
pi(x,y,z),  F(x,y;^T)  and  F(x,y;X)  over  the  cross  section  of  the  pore.  Following  the 
traditional  approach  of  porous  media  modeling, ^0,1 1  complex  compressibility  is  defined 
as 


C(z;A,t) 


_L  Pi(z) 
PO  ' 


(2.23) 


and  from  Eq.  (2.22),  with  Tqz  =  0,  is  given  by 

C(rM)  =  — 4  (y  -  <r  -  1)  F(>-t)  )  •  (2.24) 

Poc2 


The  cross-sectionally-averaged  density  can  then  be  expressed  as 


Pl(z) 


PO  C(z;Xt)  Pi(z) 


F(Xt)  -  Npr  FQ.) 

1  -  Npr 


dPi(z) 

dz 


(2.25) 


As  was  the  case  with  the  complex  density,  the  use  of  the  terminology  complex  signifies  that 
C(z;>.t)  is  a  complex  number.  The  motivation  for  averaging  the  density  over  the  cross 
section  of  the  pore  becomes  apparent  in  the  next  section. 


D.  Pressure  wave  equation  in  a  pore. 

The  continuity  equation,  Eq.  (2.7),  along  with  the  equation  of  motion,  Eq.  (2.13), 
and  the  equation  of  state  for  the  density,  Eq.  (2.25),  can  be  combined  to  yield  a  wave 
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equation  for  the  pressure  in  a  pore.  The  boundary  conditions  on  Vx  and  Vy  are  that  they 
vanish  at  the  pore  wall.  By  assuming  that  we  may  interchange  the  operation  of  z 
differentiation  and  averaging  over  x  and  y  coordinates,  the  cross  sectionally  averaged 
continuity  equation  for  the  pore  is 

-i  0)  pi(z)  +  ^  (po(z)  vz(z))  =  0  .  (2.26) 


To  make  progress  the  z  derivative  on  the  second  term  must  be  applied  to  the  right  hand  side 
of  Eq.  (2.13).  Use  of  Eq.  (2.13)  gives 

i  (po(z)  V,(Z))  =  J-  A  (^F(X)  4^) .  (2.27) 

The  derivative  of  F(X)  with  respect  to  z  occurs  because  X  depends  on  temperature  and 
hence  on  z  also.  Multiplying  Eq.  (2.26)  by  ico,  dividing  through  by  F(X),  using  the 
complex  density  of  Eq.  (2.14),  the  equation  of  state  for  density,  Eq.  (2.25),  and  Eq. 
(2.27),  the  wave  equation  for  pressure  is 

P,(z)  =  0,  (2.28) 


where 


a(A.,>.T)  = 


Tor  ^ 


-  F(Xt)/F(>,)  -  Npr  >1 

1  -  Npr 


(2.29) 


and 


C(Ut)^  =  ^P^(y-(Y-1)F(Xt))  . 


(2.30) 
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In  the  abscence  of  a  temperature  gradient,  Tqz  =  0  so  a(A,,^T)  =  0.  The  propagation 
constant  in  the  pore  is  then  given  by  ±  ^  which  is  the  usual  form  used  in  porous  media 
modeling.  lOdl  The  form  of  the  wave  equation  for  pressure  is  reminiscent  of  the  time 
analog  of  a  damped  harmonic  oscillator,  however,  here  a  and  ^  are  complex  quantities. 

E.  Specific  acoustic  impedance  and  pressure  translation  theorems. 

It  is  appropriate  to  establish  the  terminology  used  in  this  section.  The  specific 
acoustic  impedance  of  an  acoustical  medium  is  equal  to  the  ratio  of  the  total  acoustic 
pressure  and  total  particle  velocity.  For  example,  in  a  fluid  layer,  the  total  pressure  is  a 
combination  of  a  down  going  wave  and  an  upgoing  wave.  The  acoustic  impedance  is  equal 
to  the  ratio  of  the  total  acoustic  pressure  and  the  total  volume  velocity.  For  porous  media, 
the  appropriate  boundary  condition  at  a  surface  is  continuity  of  the  acoustic  impedance. 
For  adiabatic  sound,  the  characteristic  or  intrinsic  impedance  is  equal  to  po  c  where  po  is 
the  ambient  density  and  c  is  the  adiabatic  sound  speed. 

Rayleigh  developed  an  impedance  translation  theorem  for  homogeneous  fluid 
layers. The  impedance  translation  theorem  relates  the  specific  acoustic  impedance  at  one 
side  of  a  layer  to  that  at  the  other.  In  this  manner,  one  may  apply  the  theorem  as  many 
times  a  necessary  to  compute  the  specific  acoustic  impedance  at  any  surface  in  the  layered 
media. 

To  this  point,  we  have  considered  propagation  in  a  single  pore  of  infinite  length. 
Denote  by  N  the  number  of  pores  per  unit  area  in  a  porous  sample  (e.g.  Fig.  2.1b) 
consisting  of  a  parallel  combination  of  many  identical  single  pores  (e.g.  Fig.  2.1c)  of  finite 
length.  Define  by  Vzb  a  particle  velocity  for  the  bulk  porous  sample,  by  A^s  the  area  of  the 
resonator  at  the  face  ot  the  sample,  and  by  A  the  area  of  a  single  pore  in  the  sample. 
Volume  velocity  is  Ares  ^zb  =  N  A  Ares  Vz  =  Q  Ares  Vz  where  Q  =  N  A  is  the  porosity  of  the 
sample.  Thus  by  using  Vz  =  Vzb  /  in  all  of  the  single  pore  equations,  both  acoustic 
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impedance  and  specific  acoustic  impedance  are  continuous  at  each  interface  (e.g.  at  the 
stack,  heat  exchanger  interface  in  Fig.  2.1a)  between  the  porous  media. 

Several  assumptions  are  necessary  to  construct  an  impedance  translation  theorem  for 
the  stack.  First,  although  the  stack  temperature  varies  continuously  from  the  hot  end  to  the 
cold  end,  we  assume  that  the  stack  can  be  modeled  as  a  series  of  isothermal  layers.  Each 
layer  temperature  is  the  average  of  the  temperature  at  the  ends  as  indicated  in  Fig.  2.2.  The 
number  of  layers  and  layer  thickness  are  practically  unrestricted.  The  physical  parameters 
X  =  R  (po  0)  / 11)1/2  and  Xt  =  R  (po  Cp  O)  /  k)1/2  are  taken  to  be  constant  within  a  layer,  and 
are  evaluated  for  the  average  temperature  of  the  layer.  The  most  important  consequence  of 
these  assumptions  is  that  a  and  ^  in  Eq.  (2.29)  can  be  considered  constant  within  a  layer. 
The  wave  equation  for  pressure,  Eq.  (2.29),  can  then  be  solved  easily  and  an  impedance 
translation  theorem  analogous  to  Rayleigh’s  may  be  derived.  We  also  assume  that  the 
acoustic  field  equations  developed  in  Sec.  2.2A  through  Sec.  2.2D  apply  to  each  finite 
length  pore  in  the  stack. 


Z-d 


Z+d 


I 


To  +  Toz  (Z  -  d/2) 


To  +  Toz  (Z  +  d/2) 


ISOTHERMAL  SUBSECTIONS  OF  A 
THERMOACOUSTIC 
ELEMENT 


Figure  2.2  An  isothermal  subsection  of  a  thermoacoustic  element,  shown  here 
having  square  pore  capillary  tubes.  Each  subsection  is  taken  to  be 
locally  isothermal,  at  the  average  temperature  of  the  end 
temperatures.  For  the  stack  thermoacoustic  element,  an  arbitrary 
number  of  subsections  may  be  used  in  spanning  the  temperature 
difference  at  the  ends  of  the  stack.  Constant  temperature  heat 
exchangers  are  associated  with  only  one  subsection. 
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Allowing  for  these  assumptions,  the  solution  of  Eq.  (2.28)  is 
Pl(z)  =  B  expi-az)  expii  kz)  +  D  exp(-az)  exp(-i  kz). 


(2.31) 


where  B  and  D  are  constants  and 

k2  =  ;2_a2  .  (2.32) 

One  must  take  care  in  obtaining  k  since  both  ^  and  a  are  complex  quantities.  From  Eq. 
(2.13),  bulk  particle  velocity  is 

Vzb(z)  (B(ia  +  k)ej:p(-az)ejcp(ikz)  +  D(ia  -  k)exp(-az)exp(-i  kz)).  (2.33) 

topO 


From  Eq.  (2.31)  and  Eq.  (2.33),  one  may  compute  a  propagator  formalism  in  the  usual 
way. 


Pl(z-d) 


Vzb(z-d) 


=  exp (ad) 


-i  Zini5//J(kd)  VPl(z) 


co5(kd)  -  a/k5/n(kd)  -i  Zini5//j(kd) 

-i  ^  .ym(kd)  c<?^(kd)  +  ^5m(kd) 

Z-int  K 


Vzb(z) 


(2.34) 


where  Zint  =  po  w  /  (F(^)  ^2  k)  is  analogous  to  characteristic  impedance.  From  use  of  Eq. 
(2.34),  and  referring  to  Fig.  2.2,  a  right-to-left  impedance  translation  theorem  can  be 
computed 


Z(z-d)  =  Zint 


_ Z(z)  [c<?5(kd)  -  g/k  5/n(kd)]  -  i  Zjnt  5tn(kd) _ 

Zint  [co5(kd)  -I-  a/k  jm(kd)]  -  i  Z(z)  [1  +  a2/k2]  sinikd)' 


(2.35) 


and  a  left- to  right  pressure  translation  theorem  can  be  computed 


Pl(z-Kl)  =  exp(-ad)  Pi(z)  ^cor(kd)  -i- 


|§)sm(kd)) 


(2.36) 


The  impedance  translation  theorem  reduces  to  Rayleigh's  form  when  no  temperature 
gradient  is  present  on  the  stack  and  a  =  0.  Use  of  these  translation  theorems  will  be 
discussed  next. 

The  translation  theorems  in  Eq.  (2.35)  and  Eq.  (2.36)  can  be  used  to  determine  the 
impedance  and  pressure  at  all  points  in  the  resonance  tube  in  Fig.  2.1.  For  thermoacoustic 
elements  other  than  the  stack,  a  =  0  in  Eqs.  (2.35)  and  (2.36).  One  starts  with  a  known 
value  of  impedance  at  one  end  of  the  tube.  As  an  example,  the  rigid  end  of  the  tube  in  Fig. 
2.1a  a  starting  point.  Subsequent  use  of  Eq.  (2.35)  determines  Z  at  the  hot  heat  exchanger 
interface,  the  stack-hot-heat-exchanger  interface,  etc.  The  stack  may  be  broken  up  into 
several  isothermal  layers  as  indicated  in  Fig.  2.2.  Finally  the  acoustic  impedance  at  the 
acoustic  driver  can  be  computed.  For  known  acoustic  driver  response  (e.g.  perhaps 
constant  velocity  or  constant  displacement,  etc),  the  pressure  can  be  determined  at  the 
driver  location.  Then  use  us  Eq.  (2.36)  determines  the  pressure  at  all  points.  It  should  be 
noted  that  all  variables  in  Eq.  (2.35)  and  Eq.  (2.36)  are  to  be  evaluated  at  the  local  position 
of  the  thermoacoustic  element  or  subsection.  The  translation  theorem  approach  for 
computing  acoustical  quantities  is  a  superior  way  of  performing  calculations  for  resonators 
containing  many  thermoacoustic  elements.  Knowledge  of  heat  and  work  flow  is  central  to 
thermoacoustics.  In  the  next  section,  heat  and  work  flow  are  evaluated  for  arbitrary  pore 
perimeters  and  are  expressed  in  terms  of  pressure  and  specific  acoustic  impedance. 

F.  Heat  and  work  flow. 

The  time  averaged  energy  flow  H2(z)  is^ 

H2(z)  =  Q2(z)  +  W2(z)  -  Qioss(z)  (2.37) 
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where  time  averaged  heat  flow  due  to  hydrodynamic  transport  Q2(z)  is** 


Q2(z)  =  ^Re^  J(pOCpVz(x,y,z)Ti  (x,y,z)*  -  pTovz(x,y,z)Pi(z)*)  dxdy  (2.38) 


the  heat  flow  due  to  conduction  down  a  temperature  gradient  is"^ 


Qloss(z)  —  A  i^gas  Tqz  +  ( 1  —  Q)  A  Kstack  Tqz  , 


(2.39) 


and  the  time  averaged  work  flow  W2(z)  is"^ 


W2(z)  Re  ^  jvz(x,y,z)  Pi(z)*  dx  dy 


(2.40) 


Here  A  is  the  cross-sectional  area  of  the  resonance  tube  at  point  z,  Ap  is  the  cross-sectional 
area  of  a  single  pore,  is  porosity,  Kgas  and  Kstack  are  the  thermal  conductivity  of  the  gas 
and  stack,  and  Re  indicates  the  real  part  of  the  expression.  The  product  ft  A  is  cross- 
sectional  open  area  of  the  tube  at  position  z. 

To  compute  ^(z)  and  W2(z)  we  make  use  of  Eq.  (2.2 1 )  for  Ti  (x,y,z)  and  Eq.  (2.11) 
for  Vz(x,y,z)  in  Eq.  (2.38)  and  Eq.(2.39).  Resulting  expressions  are 


X  ,  ,  n  A  po  Cp  ^  rF(x,y;X)  „  ,  ^ 

(i(z)  =^— Im  - Plz(z)x 

^  J  wpo 

PiW-^  FO^y;^-  -  Np,  F(x,y;^)- 


c2  po  (3 


1  -  Npr 


-  pToW2(z) 


(2.41) 


17 


where  Piz(z)  =  ciPi(z)/dz  and  *  implies  complex  conjugation,  and 


W2(z)  = 


CIA  J_ 

2  Ap 


p^(x.y;^) 
I  wpo 


Plz(z)  Pi(z)*  dx  dy  . 


(2.42) 


Recal  the  definitions  Xj  =  Npr*^  X  and 


F(X)  =  J  F(x,y;^)  dx  dy 

P  An 


(2.43) 


The  following  general  integral  result,  which  is  proven  at  the  end  of  this  section. 


^  F(x,y;X)  F(x,y;XT)*  dxdy 


F(X)  Npr  4  F(>.t)* 

Npr  +  1 


(2.44) 


and 


l2 


^  J  F(x,y;>.)  F(x,y;X)*  dx  dy 

P  Ap 


=  Lim  I] 


Re  F{X). 


(2.45) 


are  necessary  to  evaluate  Eqs.  (2.42)  and  (2.43).  By  making  use  of  these  integrals  and  the 
thermodynamic  relation  y  —  1  =  Xq  c2  /  Cp,  heat  and  work  flow  are 


^(z)  pTo  Iml 


rPiz(z)Pi(z)*  F(Xt)*-  F^)  >, 
^  pOO)  1  +  Npr 


Q  A 
2 


IPlz(z)|2 


Im(F(XT)*  +  Npr  FjX)) 

1  -  Npr2 


(2.46) 
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are 


W2(z)  = 


Q  A 
2 


^  /Plz(z)Pl(z)* 


POCO 


(2.47) 


Equations  (2.46)  and  (2.47)  are  general  expressions  for  heat  and  work  flow  with  the 
functional  form  of  F(X)  dependent  on  the  particular  pore  geometry.  It  can  be  shown  that 
Eq.  (2.46)  and  Eq.  (2.47)  are  the  same  as  Swift's4  equations  for  the  special  case  of  parallel 
slit  geometries  and  assuming  the  Es  =  0  (appropriate  for  gas  thermodynamics)  in  Swift's 
theory. 

Since  impedance  and  pressure  translation  theorems  have  been  derived  for  arbitrary 
locations  within  the  resonator,  and  specifically  for  the  stack  element,  it  is  useful  to  express 
heat  and  work  flow  in  terms  of  the  specific  acoustic  impedance  Z(z)  and  pressure  Pi(z). 
From  the  definition  of  specific  acoustic  impedance,  Z(z)  =  Pi(z)  /(f2  Vz(z)),  and  Eq.  (2.8) 
for  V2(z), 


dPi(z) 

dz 


=  ,i°>P0PlW 
n  F(X)  z(z) 


(2.48) 


Heat  and  work  flow  are 


Q:(z)=i^PTo,P.(z)l^.™f  ‘  F(XT)VFa)-n 


a  Z(z) 


1  +  Npr 


a  A 


^  IPl(z)|2  Toz-^ 


1 


Im(F(XT)*  +  Npr  F(X)) 


0)  IF(X)QZ(z)|2 


1  -  Npr2 


.  (2.49) 
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arc 


(2.50) 

In  Eqs.  (2.49)  and  (2.50),  Pi(z)  and  Z(z)  are  global  variables  in  that  they  depend  on  the 
detailed  arrangement  of  all  elements  of  the  thermoacoustic  engine,  and  F(^)  and  F(A.t) 
depend  on  the  local  properties  of  the  gas  and  stack. 

The  remainder  of  this  section  is  devoted  to  proof  of  Eq.  (2.44).  From  Eq.  (2.12),  we 

can  write  Eq.  (2.44)  as 


W2(z)=-^  IPl(z)|2ln/  — — 
2  \nz(z) 


ii 


J_ 

Ap 


F(x,y;)t) 


+  :^V^xyF(x,y;X,T)*ldx  dy 


(2.51) 


Integration  by  parts  and  the  definition  of  F(X)  give 


Il=  F(>.)  +  ^  7^  J  ^xy  ■  (F(x,y;A,)  VxyF(x,y;XT)*)  dx  dy 

P  iArX^  An 


7^  J  VxyF(x,y;?i)  '  VxyF(x,y;XT)*  dx  dy  , 


(2.52) 


where  the  gradient  operator  in  the  x,y  plane  is  Vxy.  The  divergence  theorem  in  the  x,y 
plane  can  be  applied  to  the  second  integral  in  Eq.  (2.52).  From  the  divergence  theorem,  we 
have 


/  Vxy  ■  (F(x,y;A)  VxyF(x,y;\T)*)  dx  dy  = 

Ap 

Jn  ■  (F(x,y;X.)  VxyF(x,y;>.T)*)  dS  =  0  ,  (2.53) 
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where  dS  is  an  element  of  perimeter  on  the  pore  of  arbitrary  shape,  and  n  is  the  outward 
normal  at  the  pore  wall.  The  integral  in  Eq.  (2.53)  is  zero  by  the  boundary  condition 
F(x,y;A.)  =  0  for  X  and  y  on  S.  Thus 


Il  = 


m- 


^  J  VxyF(x,y;X)  '  VxyF(x,y;>.T)*  dx  dy 
lA-T^An 


(2.54) 


Now  we  could  have  instead  used  Eq.  (2.12)  for  F(x,y;>.)  in  Eq.  (2.44),  rather  than  for 
F(x,y;^T)*  as  we  did  to  obtain  Eq.  (2.54).  Had  we  done  this,  Eq.  (2.54)  would  be 


II  =  F(X,t)*+  ^  .J  ^xyF(x,y;X)  '  VxyF(x,y;XT)*  dxdy 


(2.55) 


Eliminating  the  common  integral  between  Eq.  (2.54)  and  Eq.  (2.55)  gives  Eq.  (2.44). 

2.3  Heat  and  Work  Flow  in  the  Short  Stack  Approximation  for  Various  Stack  Geometries 
In  this  section  heat  flow  is  compared  for  stacks  consisting  of  a  variety  of  capillary 
tubes.  In  addition  to  the  square  pore  stack  shown  in  Fig.  2.1b,  parallel  plate,  circular,  and 
equilateral  triangular  capillary  tubes  will  be  considered.  Figure  32.  shows  the  arrangment 
for  the  short  stack  approximation.  In  this  approximation,  pressure  and  specific  acoustic 
impedance  are  given  by  empty-tube  values  for  empty  tube  propagation  constant  kq, 

Z(s)  =  i  po  c  cor  (kos) ,  (2.56) 


and 


„  ,  ,  i  PO  c  Vo  co^(kos) 
— 


(2.57) 
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Use  of  Eq.  (2.56)  and  Eq.  (2.57)  in  the  heat  flow  equation,  Eq.  (2.49)  gives 


IsmHkoL)  1  +  Npr  ^ 

\  Im  {F*(?cT)  +  NprF(X))  >| 

^  Im  {F*(Xt)F(X))  (l-Npr)J 


(2.58) 


where 


P  ^  Tfe  P  Tq  c  tan  (kps) 
^0  (0  (y  -  1 ) 


In  the  inviscid  approximation  for  which  Npr  =  0,  F(A.)  =  1, 

^(z)  =  P  To  po  c  Im  F*(Xt)  (1  -  n  .  (2.60) 

Physically,  the  term  Im  F*(^x)  is  a  measure  of  the  dynamic  thermal  interaction  between  the 
gas  and  solid.  Recall  that  =  R  (po  O)  Cp  /  k)^/2  where  R  is  twice  the  ratio  of  capillary 
pore  area  to  pore  perimeter.  The  function  F(x,y;X.T)  is  a  solution  to  the  partial  differential 
equation,  Eq.  (2.12),  for  a  particular  pore  geometry  and  F(Xx)  is  the  average  of  this 
quantitiy  over  the  pore  cross-section.  According  to  Eq.  (2.60),  stacks  made  of  pores  for 
which  Im  F* (Kj)  is  a  large  value  will  result  in  the  greatest  heat  flow. 

Work  flow  is  given  generally  by  Eq.  (2.50).  No  work  is  done  in  the  region  to  the 
right  of  the  stack  in  Fig.  2.3,  which  can  be  verified  readily  by  using  the  impedance  of  Eq. 
(2.56)  in  Eq.  (2.50).  In  this  region,  pressure  and  velocity  have  standing  wave  phasing. 
Work  is  done  in  the  stack.  To  compute  the  amount  of  work  done,  we  make  use  of  the 
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impedance  translation  theorem  to  compute  the  impedance  at  the  left  side  of  the  stack.  In  the 
short  stack  approximation  we  take  kd  «  1,  hence  from  Eq.  (2.35), 


^  ^  i  PO  d  i  d  Z(s)  F(X) 

Z(s+d)  =  Z(s)  1  -  2ad  -  ^ 

1^  F(X)  Z(s)  pO  CO 


(2.61) 


(ST§  I  STACK 

V(0,t)  =  Vo  exp  -i  cot 

d 

Figure  2.3  Arrangement  for  the  short  stack  approximation. 


After  some  manipulation,  work  flow  is 


W2(z)  (Y-  l)po(odImF*iXT)co5Hkos)  fl  -  F  {F*(^T)/F;;a)}. 

2^m2(koL)  F*(;Lt)  (1  -  Npr) 


+  W2visc(z) 


(2.62) 


where  work  done  on  the  gas  due  to  viscosity  is 


^\r  /  \  •^res  j  •  7/1  X  F*(?'-) 

W2vi.sc(z)  =  L  T  ^  PO  “  ^  sinHkQs)  — - ^ 

25ZAz^(koL)  F(X)  F*(>.) 


(2.63) 


In  the  inviscid  approximation,  work  flow  is 
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W2(z)  =  (Y-  1)  PO  0)  d  Im  F*aT)  cos^kos)  (1  -  T)  . 


(2.64) 


Equation  (2.62)  and  Eq.  (2.64)  are  the  acoustic  power  absorbed  by  the  stack  and  when 
these  quantities  are  negative,  it  indicates  that  acoustic  power  is  being  produced  by  the  stack. 
Just  as  we  saw  for  heat  flow  above,  stacks  made  of  pores  for  which  Im  F*(Xt)  is  a  large 
value  will  result  in  the  greatest  work  flow. 

Work  and  heat  flow  are  to  be  compared  for  the  various  pore  geometries  shown  in 
Fig.  2.4a-d.  In  the  inviscid  short  stack  approximation,  pores  with  a  large  value  of 
Im  F(?..)*  will  have  the  greatest  heat  and  work  flows  as  indicated  by  Eq.  (2.60)  and  Eq. 
(2.64).  According  to  the  results  shown  in  Fig.  2.5,  the  parallel  slit  geometry  has  the 
largest  value  of  Im  F(X)*.  The  value  occurs  for  Xc  ~  3  2  which  allows  one  to  compute  the 
optimal  operating  frequency  from  the  relation  Xc  =  (po  co  Cp  /  k)^/2  jp  other  words,  you 
can  get  about  10%  more  heat  flow  and  work  flow  in  thermoacoustics  just  by  choosing  to 
make  your  stack  from  parallel  plates  rather  than  the  other  pore  geometries.  The  functional 
form  of  F(X)  for  the  various  pore  geometries  is  given  below. 

Swift^  and  Rott^  have  worked  out  the  parallel  slit  geometry  stack.  Parallel  slit  porous 
media  has  also  been  of  interest  in  porous  media.lld3  Denote  by  2a  the  separation  distance 
between  plates.  To  be  consistent  with  the  definition  of  the  characteristic  pore  radius  R  as 
being  equal  to  twice  the  transverse  pore  area  divided  by  the  pore  perimeter,  we  take  R  =  2a. 
The  y  axis  is  centered  at  the  mid-point  between  plates  with  the  z  axis  extending  in  the 
longitudinal  direction.  The  function  F(y;A.)  which  satisfies  the  differential  equation  for  the 
transverse  dependence  of  the  equation  of  motion,  Eq.  (2. 12),  is 


F(y;X)  =  1 


coshi'J^  X/2  y/a) 

coshi'l^  X/2) 


(2.65) 
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c)  >'  d) 


Figures  2.4  (a)  Parallel  slit,  (b)  circular,  (c)  rectangular,  and  (d)  equi¬ 
lateral  triangular  capillary  tube  geometries  considered  for  the 
short  stack  approximation.  For  example,  c)  corresponds  to 
the  stack  shown  in  Figure  lb. 


For  the  slit,  the  average  of  F(y;A,)  over  the  pore  is  defined  as  F(X)  =  l/(2a)  j  F(y;X)  dy  and 
is 


F(X)  =  1  -  -^tanhirfxrjl)) 

XVT 


(2.66) 


Rott  has  worked  out  the  cylindrical  pore*  geometry  stack.  One  type  of  porous  media 
theory  is  based  upon  propagation  in  a  single  cylindrical  cylindrical  capillary  tube.^^dl 
Denote  by  a  the  radius  of  the  cylindrical  pore.  The  characteristic  pore  radius  for  cylindrical 
pores  is  R  =  a.  The  radial  coordinate  of  a  cylindrical  coordinate  system  centered  at  the 
middle  of  the  circular  pore  is  r.  The  function  F(r;X)  which  satisfies  the  differential  equation 
for  the  transverse  dependence  of  the  equation  of  motion,  Eq.  (2.12),  is 
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Figure  2.5  Real  and  imaginary  part  of  the  F(Xt)  for  parallel  slits  (+),  circular 
pores  (closed  circles),  square  pores  (closed  squares),  and  equilateral 
triangular  pores  (closed  triangles).  Heat  flow  is  proportional  to  the 
imaginary  part  of  F*(Xt)  in  the  short  stack  approximation.  Parallel 
slit  stack  geometry  is  the  best  choice  for  optimizing  heat  flow. 


(2.67) 


FCr;X)  =  1 


Joi-4i  X  r/a) 

/o(Vrx) 


The  average  of  F(r;A.)  over  the  pore  is  defined  as  F(A.)  =  2/a2  /  F(r;X)  r  dr  and  is 


F(Ji)  =  1  -  ZzMAl 
vta,  /o(Vrx) 


(2.68) 


Previous  work  has  emphasized  thermoacoustic  stacks  consisting  of  parallel  slit  pores 
and  cylindrical  pores.  The  basic  equations  are  given  here  for  stacks  consisting  of 
rectangular  pores.  Rectangular  pores  offer  a  much  wider  range  of  applicability  than  do 
parallel  slits  or  cylindrical  pores."^  In  addition,  we  anticipate  using  square  pore  ceramic 
material  in  a  thermoacoustic  refrigerator.  Denote  by  2a  and  2b  the  length  along  the  x  axis 
and  y  axis  of  the  rectangular  pore  cross  section.  We  take  the  coordinate  system  origin  to  be 
at  the  lower  left  comer  of  the  rectangular  pore.  It  is  convenient  to  define  a  function  Ymn(^) 
as 


Ymn(>.) 


iTi^  b^m^  +  a^n^ 
(a  +  b)2 


(2.69) 


The  characteristic  transverse  dimension  equal  to  twice  the  pore  area  divided  by  the  pore 
perimeter  is 

R  =  .  (2.70) 

a  +  b 


A  series  solution  for  the  function  F(x,y;^)  which  satisfies  the  differential  equation  for  the 
transverse  dependence  of  the  equation  of  motion,  Eq.  (2.12),  is 
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F(x,y;A.) 


16  sinjmnx/ldi)  sin(nny/2h) 

7^2  mnYmn(^) 

m  n  odd 


(2.71) 


In  all  sums  given  in  this  section,  m  and  n  extend  over  positive  odd  integers.  For  a 
rectangular  pore  the  average  is  defined  by  F(X)  =  l/(4ab)  /  F(x,y;X)  dx  dy,  where  the 
integral  extends  over  uie  entire  cross  section  of  the  pore. 

F(>.)  =  ^  y  — — J -  .  (2.72) 

m2  n2  Y(X,) 

m  n  odd 


We  obtained  F(x,y;X)  in  Eq.  (2.71)  from  a  solution  by  Hanl"^  for  the  same  differential 
equation  as  Eq.  (2.12). 

The  solution  for  equilateral  triangular  pores  has  been  recently  reported  by  Stinson. 
The  pore  geometry  is  shown  in  Fig.  2.4d.  The  characteristic  dimension  is  R  =  a  /  3^(2  and 


F(X)  =  1  - 


coth 


3X^|  -i 
2 


4  i 
3X2 


(2.73) 


2.4  Conclusion 

Assume  that  one  can  compute  F(x,y;X)  from  Eq.  (2.12)  and  the  boundary  condition 
F(x,y;X)  =  0  for  x  and  y  on  the  arbitrary  pore  perimeter,  and  F(X)  by  averaging  this 
function  over  the  pore  cross-sectional  area.  All  first  order  acoustical  field  quantities  (Sec. 
2.2A-D)  and  the  second  order  energy  flux  (Sec.  2.2F)  can  be  computed  from  this  function 
and  its  companions  F(x,y;XT)  and  F(Xx).  This  general  framework  was  used  in  Sec.  2.3  to 
investigate  the  optimal  choice  of  capillary  tube  geometry  for  stacks.  In  the  inviscid,  short 
Slack  approximation,  the  parallel  slit  stack  geometry  is  best  for  maximizing  heat  and  work 
flow. 
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Impedance  and  pressure  translation  theorems  were  developed  in  Sec.  2.2E  for 
determining  these  quantities  at  all  points  in  the  resonator  shown  in  Fig.  2.1a.  With  this 
approach,  analysis  of  complicated  arrangements  of  thermoacoustic  elements  can  be 
evaluated  readily  and  in  a  unified  manner.  Work  and  heat  flows  were  expressed  in  terms  of 
specific  acoustic  impedance  and  pressure  to  take  advantage  of  these  theorems. 

Finally,  the  function  F(k)  is  also  the  key  element  of  capillary-pore  based  porous 
media  models. Factors  are  used  in  these  models  to  scale  properties  of  random 
media  to  circular  pores.  Thus  the  scaling  factors  and  methodology  of  porous  media 
modeling  can  be  readily  adapted  to  be  useful  in  thermoacoustics.  Afterall,  thermoacoustic 
elements  are  nothing  more  than  sections  of  a  porous  media,  with  the  added  complication  of 
temperature  gradients. 
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3.0  DESIGN  OF  A  THERMOACOUSTIC  HEAT  ENGINE 


The  design  of  the  thermoacoustic  engine  involves  the  development  of  criteria  for  the 
following  elements: 


1 .  Tube  length  and  operating  frequencies 

2.  Heat  exchangers 

3.  Stack 

4.  Driver 

5.  Gas  handling  system 

6.  Fluid  flow  system 

7.  Temperature  measurements 

8.  Assembly 


In  the  following  we  outline  the  design  criteria  for  tube  length  and  operating 
frequencies  and  the  design  of  the  heat  exchangers. 


3.1.  Spacing  of  plates  in  the  stack  and  the  gas  and  ambient  pressure  and  temperature 
combine  to  give  the  operating  frequency  and  resonator  tube  length  for  the 
thermoacoustic  refrigerator. 

In  Fig.  3.1,  a  is  the  semiwidth  of  each  square  of  the  ceramic  material  and 
5k  =  [2  K/(po  Cp  is  the  thermal  penetration  depth  which  is  roughly  the  distance 

heat  can  diffuse  through  the  fluid  in  time  l/o).  Here  gas  properties  are  k  =  thermal 
conductivity,  po  =  ambient  gas  density,  and  Cp  is  the  constant  pressure  heat  capacity  per 
unit  mass,  co  is  the  radian  frequency.  The  thermal  diffusivity  of  the  gas  is  K  =  k  /  po  Cp. 
The  condition  on  5k  and  2a  shown  in  Fig.  3.1  was  taken  from  Ref.  1,  pg.  1157  above  Eq. 
(60)  and  will  be  used  as  the  design  criteria.  The  operating  frequency  f  for  a  square  pore 
stack  thus  is  obtained  from  a  =  35k  or 


9  K _ 9  T] 

7C  a2  po  Cp  Npr  K  a2  po 


(3.1) 
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ACOUSTICAL  EQUIVALENCE  AMONG  SLITS  AND  SQUARES 


T" 

asl=  1.75  Sic 


Figure  3. 1  Parallel  slits  of  semiwidth  a  used  by  Swift  are  roughly  acoustically 
equivalent  to  square  of  semiwidth  1.65a.  Requirements  on 
boundary  layer  thickness  and  semiwidth  are  shown. 


Gas  density,  pressure,  and  temperature  are  related  by  the  ideal  gas  relation 


P  = 


_ P_ 

T  R/|irnw 


(3.2) 


where  R  =  8.3143  joules  mole"'  Kelvin'^  is  the  gas  constant  and  Pmw  i  s  the  average 
molecular  weight  in  kilograms  /  mole  and  for  Helium  (Argon)  jimw  =  4  /  1000  (40  /  1000) 
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kilograms  /  mole.  For  a  monatomic  gas  Cp  =  (5/2)  R  /  |imw  From  Ref.  2,  pg  230-231, 
viscosity  of  Helium  and  other  monatonaic  gases  follows 


ri  =Ti’  (TyT')s  ,  (3.3) 

where  for  Helium  (Argon)  we  may  take  T)'  =  1.887  *  lO'^  (2.117  *  lO'^)  Pa  sec  for  T'  = 
273.15  Kelvin  and  s  =  0.6567  (0.81 1)  for  the  temperature  range  33.15  Kelvin  to  1073.15 
Kelvin  (293  to  373  Kelvin).  For  a  monatomic  gas  the  Prandtl  number  is  (see  Susalla's 
thesis.  Ref.  3,  and  references  therein)  Npr  =  T]  Cp  /  K  =  2/3  so  we  may  calculate  the  thermal 
conductivity  K  as 

K  =  T)  Cp  /  Npr  =  T1  (5/2)  R  /  4mw  /  (2/3)  =  K’  (TyT)s  ,  (3.4) 

where  for  Helium  (Argon)  k'  =  0.1471  (0.01650)  Watts  /  (meter  Kelvin)  at  T  =  273.15 
Kelvin  and  s  =  0.6567  (0.81 1)  for  the  temperature  range  33.15  Kelvin  to  1073.15  Kelvin 
(293  to  373  Kelvin).  The  sound  speed  for  an  ideal  gas  is  given  for  reference  here  as 

c  =  =c'(T7r)0.5  (3.5) 

1  H-m  w 


where  for  Helium  (Argon)  c'  =  972.8  (307.6)  m/s  at  T  =  273.15  Kelvin. 

Consideration  of  the  operating  frequency  for  Helium  (Argon)  gas  we  have  from  Eqs. 
(3. 1-3.3) 


fHe  (fAr)  - 


TqI.6567 
236  a2  po 


TqI-SH 
5005  a^  PO 


(3.6) 


Lets  work  at  an  ambient  temperature  of  Tq  =  300  K  and  an  ambient  pressure  po  =  3  atm  =  3 
atm  *  101325  Pa  /  atm  =  303975  Pa,  and  use  the  square  pore  media  having  a  =  0.77  mm  so 
that  we  get  an  operating  frequency  of  fHe  =  300  Hz  and  fAr  =  34  Hz.  The  operating 
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wavelength  is  Xhc  =  c  /  fne  =  3.4  m  and  Xat  =  c  /  Fat  =  9.5.  For  a  pressure  maximum  at 
both  ends  of  the  tube,  Ltube  =  ^He  /  2  =  1.7  m  for  Helium  and  Lmbe  =  ^Ar  /  2  =  4.75  m  for 
Argon.  Tables  3.1  and  3.2  summarizes  these  results.  I  recommend  using  Helium  and  the 
criteria  given  in  the  table  for  Helium. 


TABLE  3. 1  Design  Criteria  for  Helium 


QUANTITY 

VALUE 

CRITERIA 

a  =  semiwidth  of  pore 

0.77  mm 

High  porosity,  several  lengths,  uniform 

PO  =  ambient  pressure 

3  atm 

nominal  value  =  303975  Pa 

To  =  ambient  temp. 

300  K 

nominal  value 

po  =  ambient  density 

0.488  kg  /  m3 

nominal  value 

fHe  =  operating  freq 

300  Hz 

Boundary  layer  thickness  =  a  /  3 

c  =  adia.  sound  speed 

1020  m/s 

Ambient  temperature 

^Hc  =  wavelength 

3.4  m=  11'2" 

fHe  and  c 

Ltube  ~  ^He  /  2 

1.7  m  =  57" 

operating  wavelength,  velocity  node  at  each  end 

TABLE  3.2  Design  Criteria  for  Argon 

QUANTITY 

VALUE 

CRITERIA 

a  =  semiwidth  of  pore 

0.77  mm 

High  porosity,  several  lengths,  uniform 

PO  =  ambient  pressure 

3  atm 

nominal  value  =  303975  Pa 

To  =  ambient  temp. 

300  K 

nominal  value 

Po  =  ambient  density 

4.88  kg  /  m3 

nominal  value 

fAr  =  operating  freq 

34  Hz 

Boundary  layer  thickness  =  a  /  3 

c  =  adia.  sound  speed 

322.4  m/s 

Ambient  temp>erature 

^Ar  =  wavelength 

9.5  m 

fAr  and  c 

Ltube  =  ^Ar  /  2 

4.75  m  =  157" 

operating  wavelength,  velocity  node  at  each  end 
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3.2.  Heat  Exchangers. 

Heat  carried  by  the  stack  is  deposited  at  the  ends  where  it  is  picked  up  by  the  heat 
exchangers.  Design  criteria  for  heat  exchangers  are  as  follows: 

1 .  Height  H  is  determined  by  the  peak  acoustic  displacement  in  the  gas.  The  acoustic 
displacement  amplitude  is  roughly  ^(z)  =  v^fz)  /  -io)  =  pi(z)  /  ( -ico  p  c)  and  hence  depends 
on  the  driving  amplitude.  Here  p  and  c  are  complex  quantities  as  a  result  of  dissipation  in 
the  stack  and  velocity  should  be  multiplied  by  q/Q  where  q  is  tortuosity  and  Q  is  porosity. 
For  a  rough  estimate  of  ^  we  proceed  by  supposing  pi(z)  =  0.05  po  =  15200  Pa,  po  = 
0.488  Kg/m^,  c  =  1020  m/s,  and  an  operating  frequency  of  300  Hz  to  get  ^  =  1.6  cm. 

2.  The  hot  heat  exchanger  (HHE)  is  located  near  a  velocity  node  and  the  cold  heat 
exchanger  (CHE)  is  closer  to  a  velocity  antinode.  As  commented  by  Swift  in  Ref.  1,  pg. 
1168,  the  HHE  must  pick  up  more  heat  than  the  cold  and  hence  it  should  be  larger. 
Because  the  HHE  is  close  to  a  velocity  node  viscous  dissipation  will  be  less  than  the  CHE 
so  the  HHE  can  have  a  larger  H  and  smaller  Ta  than  the  CHE. 

Combining  1.  and  2. 1  estimate  Hhhe  =  2.0  cm  and  Hche  =  1-5  cm.  For  the  HHE, 
Ta  =  0.5  mm  and  Tc  =  0.25  mm.  For  the  CHE,  Ta  =  0.7  mm  corresponding  to  the  plate 
spacing  in  the  stack  and  Tc  =  0.25  mm.  In  both  the  CHE  and  HHE,  Tq  is  the  approximate 
stack  wall  thickness  also.  Referring  to  Fig.  3.2  the  number  of  each  aluminum  and  copper 
plates  needed  is 


nplates  =  160  mm  /  (Ta  +  Tc)  (3.7) 

where  Ts  are  the  thickness  of  the  aluminum  and  copper  plates.  For  the  CHE  this  gives 
Opiates  =  120  and  for  the  HHE  Opiates  =  215.  If  we  close  every  10  channel  to  flow  fluid 
through  as  indicated  by  the  closed  channels  in  Fig.  3.2,  the  number  of  open  plates  is  Oopen 
=  nplates  *  0-9.  Hence  the  porosity  is 

Q  =  Ta  *  0.9  Opiates  /  1 60  mm  (3.8) 
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Figure  3.2.  Heat  exchangers  before  electroplating. 
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and  gives  ficHE  =  0-67  and  Ohhe  =  0-60.  By  comparison  the  stack  porosity  is  fiisTACK  = 
0.73.  Table  3.3  summarizes  the  geometry  of  the  HE's.  Table  3.4  summarizes  the  number 
of  2'X4'  sheets  of  aluminum  and  copper  needed.  Heat  exchangers  are  made  by  the 
following  steps: 

1 .  Obtain  aluminum  and  copper  plates  of  the  dimensions  given  in  Table  2. 1 .  Epoxy 
these  plates  together  as  indicated  in  Fig.  2.  The  aluminum  plate  at  the  tenth  slot 
will  have  a  height  of  Hal  =  H  -  16Tc  =  H  -  4  mm  so  that  a  channel  can  be  formed 
through  which  fluid  or  heating  wires  may  run. 

2.  Machine  the  composite  to  the  form  shown  in  Fig.  3.3  after  the  epoxy  has  dried. 

3 .  Electroplate  the  machined  exchangers  with  copper. 

4.  Machine  the  top  and  bottom  of  the  exchangers  until  the  aluminum  just  shows  in 
the  open  channels.  Do  not  machine  open  the  closed  channels. 

5.  Observing  the  locations  of  the  closed  channels,  go  to  the  side  and  machine 
through  the  electroplate  to  expose  the  aluminum.  There  will  be  2  mm  on  top  and 
bottom  of  the  closed  channels  that  are  electroplate  so  be  careful  not  to  remove  it. 

6.  The  exchanger  is  ready  to  be  dumped  into  an  acid  bath  to  remove  the  aluminum. 
Presumably  the  acid  bath  will  put  the  aluminum  into  solution,  leaving  the  copper 
relatively  unharmed. 


TABLE  3.3  Heat  exchanger  geometries  in  Metric  and  English  Units. 


Heat  Exch. 

H 

Hal 

Ta 

Tc 

^alshorts 

a 

(mm) 

(mm) 

(mm) 

(mm) 

■ 

COLD  CHE 

15 

11 

0.7 

0.25 

170 

17 

0.67 

HOTHHE 

20 

16 

0.5 

0.25 

215 

22 

m 

Heat  Exch. 

H 

Hal 

Ta 

Tc 

IRH 

^alshorts 

a 

(lO^in) 

(lO^in) 

m 

(mils) 

(mils) 

(mils) 

(mils) 

■ 

COLD  CHE 

591 

433 

28 

10 

170 

17 

0.67 

HOTHHE 

787 

630 

20 

10 

215 

22 

■ 
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TABLE  3.4  Number  of  2'  by  4'  sheets  of  A1  and  Cu  needed  for  3  sets  of  HE's. 


Heat  Exch. 

#20  mil  Ai 

#28  mil  Al 

#  10  mil  Cu 

2'  X  4'  sheets 

2'  X  4'  sheets 

2'  X  4'  sheets 

CHE 

3 

3 

HHE 

3 

3 

TOTALS 

3 

3 

6 

TOP  VffiW  AFTER  MACHINING  AND  ELECTROPLATING: 

SOME  CHANNELS  ARE  CLOSED  TO  ALLOW  FOR  FLUID  FLOW 
OR  FOR  INSERTION  OF  HEATING  WIRES 

Figure  3.3.  Top  view  of  heat  exchangers  after  electroplating  and  machining. 
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Abstract 


The  frequency  response  of  a  ihermoacoustic  prime  mover  (consisnng  of  a  parallel  plate  prime 
m.over  stack  and  heat  exchangers  housed  in  a  resonant  tube)  has  been  measured  as  a  funcrion  of  the 
m^ean  gas  pressure  and  temperature  gradient  across  the  prime  m.over  stack.  The  frequency 
response  of  rvo  simpler  resonator  configurations  was  also  mjeasured  in  the  absence  of  an  applied 
gradient.  The  quality  factor  Q  and  resonar, ce  frequency  can  be  determined  from  thio  response.  As 
the  temperaruie  gradient  is  Lncreased,  the  Q  increases  indicating  a  decrease  in  attenuation  across  the 
stack.  At  sufficiently  large  temiperature  differences  (-  300  K),  the  resonator  goes  into  self- 
oscillation  Lndicating  negative  attenuation.  Measuremenu  are  reponed  for  tiie  prime  mover  and  the 
siujpler  resonators  in  the  absence  of  a  temiperature  gradient  for  helium  and  argon  at  pressures 
ranging  from  170  -  500  .kPa.  Measuiem.ents  are  also  leponed  for  the  prime  miover  in  the  presence 
of  tem.perature  gradients,  ranging  &om  zero  to  that  required  for  onset  of  self-oscillation,  for  helium 
at  the  sarnie  pressures.  The  results  are  explained  in  terms  of  a  counterpropagating,  plane  wave 
analysis,  based  on  techniques  commonly  used  in  porous  media  investigations.  In  general,  the 
agreemient  between  theory  and  experiment  is  good,  especially  for  the  simpler  resonator 
conrigurations.  For  the  prime  mover  configuration,  there  is  considerable  disagreemient  between  the 
predicted  and  measured  Q  at  low  temperature  gradients. 


INTRODUCTION 


Swift^  and  others^’^  have  developed  theories  that  can  be  used  to  predict  the  onset  of 
oscillation  and  the  quality  factor  Q  of  a  th_rmoacoustic  prime  mover.  Their  models  are  based  on  a 
thermodynamic  approach  that  considers  energy  transfer  and  dissipation  on  the  exposed  surfaces. 
This  treatment  agrees  with  experimental  results  and  includes  all  the  applicable  physics.  It  does 
have  one  disadvantage.  As  formulated,  it  is  not  simple  to  transfer  extensive  research  of  the 
acoustic  properties  of  porous  materials  to  performance  predictions  of  thermoacoustic  devices.  One 
goal  of  this  article  is  to  develop  such  a  formalism,  which  will  be  especially  attractive  when 
considering  the  use  of  alternate  geometries  in  the  prime  mover  stack. 

In  the  following,  we  have  attempted  to  apply  the  standard  treatment  of  acoustic  waves  in 
porous  media  and  propagation  in  a  tube  with  multiple  boundaries  to  predictions  of  prime  mover 
performance.  Except  for  one  equation  for  heat  flow  in  the  presence  of  a  temperamre  gradient,  all 
equations  and  their  description  can  be  found  in  standard  acoustics  textbooks."^  As  one  might 
expect,  all  the  needed  expressions  without  a  temperature  gradient  can  be  found  in  the  works  of 
Lord  Rayleigh.^ 

1.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

We  measured  the  frequency  response  of  a  thermoacoustic  prime  mover  as  a  function  of  the 
mean  gas  pressure  and  temperature  gradient  across  the  prime  mover  stack.  The  quality  factor  Q 
and  resonance  frequency  can  be  determined  from  the  response.  The  temperature  gradients  ranged 
from  zero  to  that  required  for  the  onset  of  self  oscillation.  To  gain  confidence  in  our  ability  to 
measure  the  frequency  response,  we  first  measured  the  response  of  two  simpler  resonator 
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configurations:  an  empty,  rigidly  terminated  resonant  tube,  and  the  empty  tube  plus  the  ambient 
heat  exchanger  These  configurations,  and  that  of  the  prime  mover,  are  discussed  in  Section  A. 
Temperature  control  and  measurement  are  discussed  in  Section  B,  followed  by  a  discussion  of  the 
data  acquisition  procedure  in  Section  C. 

A.  Resonator  Configurations 

We  measured  the  frequency  response  of  three  different  resonator  configurations:  an  empty, 
rigidly  terminated  resonant  rube;  the  empty  tube  plus  the  ambient  heat  exchanger,  and  the  prime 
mover.  The  empty  resonator  is  made  from  two  3.82  cm  ID  (inner  diameter)  copper  tubes, 
separated  by  a  2.19  cm  long  brass  section  having  a  slightly  larger  ID.  This  brass  section  is  called 
the  ambient  heat  exchanger  container.  The  lengths  of  the  two  sections  of  copper  tube  are  1 1.43  cm 
and  87.95  cm.  One  end  of  each  section  is  fitted  with  a  flange  that  auows  them  to  be  soldered  to  the 
ambient  heat  exchanger  container.  A  0.64  cm  thick  copper  cap  is  soldered  to  the  other  end  of  each 
section,  forming  a  closed,  rigid,  termination.  The  cap  on  the  longer  section  houses  a  1.9  cm 
diameter  electret  driver  and  a  0.59  cm  diameter  electret  microphone.  The  driver  and  microphone 
are  flush  mounted  and  sealed  in  the  end  cap  with  epoxy.  The  end  cap  is  epoxied  to  the  section  of 
tube.  Even  vrith  the  driver  and  microphone  mounted  in  the  end  cap,  our  measurements  show  that  it 
still  behaves  as  a  rigid  termination.  The  total  length  of  the  resonator,  including  the  heat  exchanger 
container,  is  101.56  cm  (internal  dimensions).  The  resonator  is  connected  to  a  gas  handling 
system,  through  a  fill  tube  in  the  ambient  heat  exchanger  container.  This  connection  allows  for 
evacuation  of  the  resonator  before  pressurization  with  gas  (helium  or  argon)  and  a  system  venL  A 
dial  pressure  gauge  and  an  OMEGA  Model  PX304-150AV  pressure  transducer  are  connected  to 
the  fill  line  to  sense  the  mean  pressure  inside  the  tube.  The  driver  input  signal  is  provided  by  an 
HP  3325A  function  generator  and  amplified  by  a  Techron  Model  7520  power  amplifier.  The 
microphone  output  signal  is  amplified  by  a  preamplifier  with  a  gain  of  100. 

The  second  configuration  differs  from  the  first  only  in  that  the  ambient  heat  exchanger  is 
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mounted  In  the  ambient  heat  exchanger  container.  The  ambient  heat  exchanger  consists  of  two 
stacks  of  parallel  plates.  Each  stack  is  composed  of  25,  0.C45  cm  thick,  1.02  cm  long  copper 
plates.  The  gap  benveen  each  plate  is  0.104  cm.  The  two  stacks  are  mounted  in  a  3.82  cm  BD 
copper  tube  such  that  there  is  a  0.15  cm  gap  between  them.  This  double  stack  heat  exchanger  is 
not  necessary,  a  single  stack  would  suffice  for  these  measurem^ents.  But,  the  double  stack  will  be 
needed  in  a  future  phase  of  this  research. 

The  third  configuration,  the  prime  m.over,  differs  from  the  second  by  replacing  the  shoner 
section  of  copper  tube  used  in  the  first  two  configurations  with  a  different  section  This  section 
consists  of  a  nickel  heater  section,  a  nickel  hot  heat  exchanger,  and  a  stainless  steel  prime  mover. 
It  is  designed  to  withstand  the  high  temperatures  required  to  reach  onset  of  self  oscillation.  The 
heater  section  consists  of  a  3.82  era  ID,  5.762  cm  lon<^.  nickel  tube.  One  end  of  the  tube  is  rigidly 
capped  and  drilled  to  accommodate  a  thermocouple  probe,  used  to  sense  the  hot  heat  exchanger 
temperature.  The  hot  heat  exchanger  is  mounted  at  the  other  end  of  the  heater  section.  This  heat 
exchanger  consists  of  25,  0.045  cm  thick,  0.762  cm  long  nickel  plates.  The  gap  between  each 
plate  is  0.104  cm.  Between  each  plate  there  are  304  stainless  steel  spacers,  0.102  cm  long  and 
0.031  cm  in  diameter.  The  distance  from  the  closed  end  of  the  heater  section  to  the  beginning  of 
the  hot  heat  exchanger  is  5.58  cm.  The  prime  mover  stack  consists  of  35,  0.28  cm  thick  304 
stainless  steel  plates,  3.50  cm  long  and  spaced  by  0.077  cm.  This  stack  is  housed  within  a  thin 
walled  (0.05  cm  thick)  stainless  steel  tube.  One  end  of  this  stainless  steel  tube  is  welded  to  the 
nickel  tube  such  that  the  hot  heat  exchanger  and  prime  mover  stack  are  in  contact  The  other  end  of 
the  stainless  steel  rube  is  fitted  with  a  flange  so  that  the  entire  section  can  be  soldered  to  the  ambient 
heat  exchanger  container.  Summarizing  the  prime  mover  configuration,  it  consists  of  a  5.58  cm 
long  nickel  tube,  a  0.762  cm  long  hot  heat  exchanger,  a  3.42  cm  long  prime  mover  stack,  a  2.18 
cm  long  double  stack  ambient  heat  exchanger,  and  an  87.95  cm  long  copper  tube.  The  total  lenguh 
of  the  prime  mover  is  99.9  cm.  The  configuration  is  shown  in  Fig.  1. 


B.  Temperature  Control  and  Measurement 

The  nature  of  the  experiment  requires  that  a  temperature  difference  be  established  across  the 
prime  mover  stack  and  maintained  precisely  (within  1  K).  It  is  also  necessary  to  eliminate 
temperature  gradients  elsewhere  in  the  resonator.  Temperature  control  of  the  prime  mover  is 
achieved  as  follows.  An  Omega  Engineering  Model  HBA  202040  nozzle  heater  surrounds  the 
nickel  heater  section,  or  hot  end.  Electrical  power  is  supplied  to  the  heater  through  a  variable  AC 
transformer.  The  hot  end  with  its  heater  is  surrounded  by  insulation  to  reduce  heal  loss  to  the 
room  and  to  help  maintain  a  uniform  temperature  distribution  along  it  The  ambient  heat  exchanger 
is  maintained  at  near  room  temperature  by  circulating  water  through  a  jacket  which  surrounds  it. 
The  temperature  of  the  water  is  maintained  with  a  Neslab  Model  RET- 110  constant  temperature 
bath.  After  the  water  leaves  the  jacket  it  circulates  through  flexible  plastic  tubing  that  is  wrapped 
around  the  majority  of  the  length  of  the  longer  section  of  copper  tube.  This  section  of  rube  is  also 
surrounded  by  styrofoam  insulation  to  help  maintain  a  uniform  temperature  distribution  along  its 
length. 

The  temperature  of  the  prime  mover  is  monitored  by  four  thermocouples.  A  type  K 
thermocouple  is  in  contact  with  the  center  of  the  hot  heat  exchanger  to  sense  the  temperature  of  the 
hot  end.  Three  type  E  thermocouples  were  glued  to  the  top,  middle,  and  bottom  of  the  cold  end  to 
sense  the  temperature  along  that  section.  The  reference  temperature  for  the  whole  system  is  found 
by  using  a  4-wire  resistance  measurement  of  a  thermistor  mounted  on  an  aluminum  isothermal 
block.  For  the  first  two  configurations,  the  hot  end  is  replaced  by  the  shorter  section  of  copper 
tube  and  the  type  K  thermocouple  is  not  used. 

C.  Data  Acquisition 

The  quality  factor  of  the  various  resonator  configurations  is  determined  by  measuring  the 
steady  state  frequency  response.  This  measurement  is  accomplished  by  driving  the  resonator  at 
frequencies  near  resonance  and  measuring  the  steady  state  amplitude  of  the  microphone  output 
signal  with  a  Stanford  Research  Model  SR-530  lock-in  amplifier.  The  Q  is  determined  by 
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performing  a  least  squares  fit  of  the  data  to  the  ideal  response.  Data  acquisition  is  performed  by  a 
Standard  286  penonal  computer.  Referring  to  Fig.  2,  the  computer  communicates  with  the  lock-in 
amplifier,  the  HP  3457A  digital  multimeter  and  the  HP  3325A  function  generator  through  a  GPIB 
interface.  Through  the  execution  of  the  controlling  program,  a  source  signal  is  supplied  by  the 
function  generator  to  the  electnet  driver.  The  output  voltage  from  the  microphone  is  amplified  by  a 
preamplifier  and  sent  to  the  lock-in.  The  output  of  the  lock-in,  as  well  as  all  other  data  signals  of 
interest,  are  fed  to  an  HP  3457A  multimeter.  The  output  of  the  multimeter  is  recorded  by  the 
computer.  A  Kikusui  Model  COS6100A  oscilloscope  is  used  to  monitor  the  driver  input  signal 
and  the  output  signal  from  the  preamplifier. 

Before  data  acquisition  is  started,  the  resonator  is  evacuated  and  filled  with  helium  or  argon 
to  the  desired  pressure.  When  data  acquisition  begins,  the  program  records  the  mean  gas  pressure 
within  the  resonator.  Next,  the  approximate  resonance  frequency  and  the  half  power  bandwidth 
are  entered  into  the  computer,  which  then  determines  Lhe  start  and  stop  frequencies  and  the 
frequency  increment.  The  program  sets  the  driving  frequency  and  measures  temperatures, 
frequency,  and  the  output  of  the  lock-in.  The  program  then  increments  the  frequency  and  repeats 
the  process.  The  time  required  to  measure  the  frequency  response  is  approximately  5  minutes. 

The  temperature  along  the  length  of  the  resonator  is  held  uniform  during  the  measurements 
on  the  first  two  configurations.  For  the  prime  mover  configuration,  the  temperature  of  the  hot  end 
is  set  by  adjusting  the  electrical  power  supplied  to  the  heater.  After  the  hot  end  temperature  has 
reached  thermal  equilibrium,  data  acquisition  is  initiated.  The  magnitude  of  the  input  electrical 
power  is  determined  by  taking  the  product  of  the  voltage  and  current.  The  current  is  measured 
through  a  Tektronix  Model  P6021  current  probe.  The  heater  impedance  is  almost  purely  resistive, 
so  there  is  no  need  to  measure  the  phase  angle  between  the  heater  voltage  and  currenL 


II.  RESULTS  AND  ANALYSIS 
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Figure  2:  Schematic  diagram  of  the  instrumentation  used  in  the  experiment. 


the  boundary  conditions.  We  analyze  the  case  of  no  temperature  gradient  in  Section  A  below. 

Except  in  regions  where  plates  are  present,  the  propagation  constants  are  computed  from  Shield's 

work.^  In  the  regions  where  plates  are  present,  the  propagation  constants  are  computed  from  the 

7  9 

theory  of  sound  propagation  through  porous  media.  '  Pertinent  results  from  porous  media  theory 
are  discussed  in  Appendix  B.  Those  unfamiliar  with  this  t>'pe  of  analysis  are  advised  to  read 
Appendix  B  before  proceeding.  The  case  where  temperature  gradients  are  present  is  analyzed  in 
Section  B,  w’here  we  modify  the  analysis  to  include  the  effects  of  a  temperature  gradient. 

A.  Analysis  of  Data  Taken  with  No  Temperature  Gradient 

In  the  empty  resonator  case,  we  set  k2,  Z2^  k^,  Z3,  k4,  Z4,  k5,  Z5  equal  to  k^,  z^.  k^  is 
computed  from  Shield's  work^  as  k^  =  co/ci  -  ia^.  The  impedance  zj  is  computed  as  z^  =  Cdpo/k^^. 
The  length  of  the  empty  resonator  is  1.0156  m.  The  measured  and  predicted  Q's  and  resonance 
frequencies  are  shown  in  Table  I  for  the  resonator  filled  with  argon  and  helium  at  various 
pressures.  The  experimental  and  calculated  results  are  in  agreement  to  within  the  program's  ability 
to  extract  fg  (±0.5  Hz)  and  Q  given  the  uncertainty  in  temperature. 

For  the  resonator  plus  the  ambient  heat  exchanger,  Ij  is  set  at  87.95  cm;  I2  at  90.13  cm.  The 
heat  exchanger  consists  of  parallel  plates  with  a  gap  of  1.02  mm.  The  porosity  is  0.667. 
Attenborough's  theory  is  used  to  compute  Z2  and  k2.  k3  -  k5  and  Z3  -  Z5  are  set  equal  to  k^  and  zj, 
respectively.  The  tube  length  is  1.0156  m.  Results  are  shown  in  Table  H.  Again  experimental 
and  theoretical  results  are  in  close  agreement.  Note  that  tlic  sign  on  the  error  in  Q  has  changed 
signs  for  helium  suggesting  that  more  viscous  attenuation  is  predicted  than  actually  observed. 

Finally,  for  the  case  of  the  prime  mover  with  no  temperature  gradient,  I3  is  set  at  93.635  cm, 
I4  at  94.397.  The  full  length  of  the  tube  is  99.9  cm.  The  stack  consists  of  parallel  plates  with  a 
gap  of  0.79  mm  and  a  porosity  of  0.76.  This  porosity  was  multiplied  by  the  factor  of  open  stack 
area  that  is  blocked  by  the  heat  exchanger  (0.7)  to  give  an  effective  porosity  of  0.53. 
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Data  was  collected  for  the  empty  resonator,  the  resonator  wth  one  heat  exchanger,  and  the 
resonator  plus  prime  mover  stack  and  two  heat  exchangers.  The  Q  of  the  resonator  was 
determined  by  fitting  the  steady  state  frequency  response  to  a  standard  resonance  equation 

.2] 


A  (a)  = 


t1/2 


1 


(1) 


where  co  is  the  angular  frequency  of  the  drive  and  cuo  the  resonance  angular  frequency 
(frequency  of  maximum  amplinide). 

To  model  the  experiment,  a  closed  tube  with  counterpropagating  plane  waves  is  assumed  as 
shown  in  Fig.  3.  In  each  region,  the  acoustic  pressure  of  the  incident  and  reflected  waves  is 
expressed  as  p-”^  =  and  p^^^  =  respectively.  Continuity  of  pressure  and 

volume  velocity  is  assumed  at  each  interface.  Continuity  of  volume  velocity  is  written  in  the  form 
(P-(i)-Pj(i))/zj  =  (Pj(2)-Pj(2))/z2,  where  z  is  called  the  impedance  of  the  region  and  defined  in 
Appendix  A.  At  the  driving  end  (x=0),  u(0,  t)  =  ico/e'“^  where  /  is  the  effective  driver 
displacement  (assumed  to  be  independent  of  load).  The  total  complex  acoustic  pressure  at  the 
driving  end  Pjotai  is  Pi^'^  +  P/^^.  The  amplitude  of  the  sound  at  the  driving  end,  where  the 
response  is  measured,  is  the  absolute  value  of  Pto^ai-  The  expression  for  Pjotal  derived  in 
Appendix  A.  As  the  frequency  assumed  for  the  calculation  is  incremented,  Ptoul  through  a 
maximum  giving  the  resonance  frequency  fg.  Two  methods  were  used  to  find  the  values  of  Q. 
From  Eq.  (1),  values  of  A  computed  for  any  two  frequencies  co  can  be  used  to  solve  for  Q.  Values 
of  Q  were  computed  from  the  amplitudes  at  adjacent  frequency  steps  for  all  frequencies  near  fg 
using  Eq.  (1)  and  then  averaged  The  second  method  is  to  compute  values  of  A  for  frequencies 
spanning  resonance  and  then  perform  a  least  squares  fit  to  Eq.  (1)  using  Q  as  a  one  of  the 
parameters. 

To  calculate  P[ot2i»  ^  ^  specified  are  the  propagation  constants  in  each  region  and 


57 


58 


Figure  3:  Diagram  showing  the  subdivision  of  the  prime  mover  used  for  the 
counterpropagating  plane  wave  ajialysis. 


Attenborough's  theory  is  again  used  to  compute  k3  and  Z3.  Results  are  shown  in  Table  in. 

B.  Analysis  of  Data  for  a  Prime  Mover  Below  Onset 

The  expressions  for  the  propagation  constant  and  impedance  derived  in  Appendix  B  must  be 
altered  when  temperature  gradients  are  present.  To  find  these  new  expressions,  we  take  the  same 
approach  as  in  Appendix  B.  The  only  difference  being  that  we  use  a  theoretical  approach  presented 
by  Swift^  to  solve  for  the  acoustic  velocity  and  temperature  distribution  across  a  slit  in  the  presence 
of  a  temperature  gradient.  Using  these  results,  we  find  expressions  for  the  complex  density  and 
compressibility  and  finally  k  and  z.  Once  we  have  k  and  z  we  can  find  Ptoial- 

The  complex  density  is  the  same  as  that  derived  for  the  no  temperature  gradient  case  and 
derived  in  Appendix  B.  The  result  is 


Pc  =  pf  1- 


tanh[(l  +  i)a/Sv]y^ 

(1  +  i)a/Sv  I 


Swift  writes  for  the  appropriate  heat  balance  equation, 


(2) 


ico0  +  Ui 


dx  ) 


-  =  K 


d^Q 


(3) 


where  pf  is  the  density  of  the  gas;  Cp  the  isobaric  specific  heat  capacity  (per  unit  mass),  9  the 
variation  in  temperature  due  to  the  acoustic  wave,  uj  the  acoustic  velocity  in  Lhe  direction  of 
propagation,  dT^/dx  the  temperature  gradient  across  the  stack,  T^  the  mean  temperature,  .8  the 
coefficient  of  volume  expansion,  p  the  acoustic  pressure,  and  Xf,  the  coefficient  of  thermal 
conduction.  Again  from  Swift, 
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i  dp 

Ui  =-i— 

copf  dx 


1- 


cosh  ^(l+i)y/5vj 


cosh  ( (l+i)a/5v 


(4) 


Substituting  this  expression  for  into  Eq.  (3)  gives 


Ptcp  jicae  4-i-  [  1-  c°4a^i)y/5vl|  d^l . 

\  copf  dx  ^  cosh[(l+i)a/5v]/ 


ico  TmPp  =  Xh 


8^9 

3y2’ 


(5) 


where  5^  =  (2co/v)^/2;  v  is  the  kinematic  viscosity;  and  a  is  the  semi-width  of  the  slit.  Swdft 


identifies  a  solution 


0=Mp._l 


\  CTcosh[(l+i)y/5v]  ^dpdXm 

i-  *  ~r  T I  «  « 


PfCp  PfCO^  \  (a-l)cosh[(l+i)a/5v]/ 

TmP  ,  (dp/dxXdTn/dx)  ^  Eg  fy  j  cosh  [(1+i)  y/5j 
PfCp  (cr-l)pj(^2  I  fk  /  cosh  [(1+i)  a/Sj 


(6) 


where  a  =  v/k^  is  the  Prandtl  number,  and  the  other  new  terms  are  defined  in  Table  IV.  As 
explained  in  Appendix  B,  the  condensation  of  the  gas  in  the  slit  can  be  expressed  as 

s=P/Pni-6^m- 


Substituting  Eq.  (6)  into  Eq.  (7),  factoring  out  p/7Pm>  setting  =  1,  and  recognizing  that 
Pn/pfCpTj^  =  (Y-  D/yfor  an  ideal  gas,  gives 


_  P 
YPm 


S  =  ^-(1  -i-(y-  1) 


Cp  8p/9x  8T, 


m 


ci>2  p  ax 


1  +  1 


1 - ?? 

(a-DfJ 


|cosh[(l  +i)y/5j 


(a  - 1)  fj  cosh[(l  +  i)a/sjj 


+  i,C0ih[a±iW5^l  (8) 

cosh!.(l  -I-  i)a/5j j 


Integrating  the  condensation  a  CI"0  j  s  uhe  slit  gives 
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_p_  Cp  3p/3x  9Tn,  ^  £,  fy  jtanh[(l+i)a/5j  ^  j^^tanh[(l  +  i)3/5j\^  ,g. 

^  YPni\  ^  "  aj2  P  dx[  ^  (cr-l)fKi  (l+i)a/5>:  J  (l+i)a/5K  /' 

We  use  ZwiJcker  and  Kosten's  ^  definition  of  complex  compressibility  (C(a))  =  ?/p  =  s/p)  to  get 

/Pm\  co2  P  ax  [  I  (ct-l)fK/  (l+Oa/Sx  J  (l+i)a/5K  / 

For  plane  traveling  waves,  3p/9x  =  ±  ilc(Tj^)p.  Assuming  that  the  imaginary  part  of  the 
propagation  constant  is  small  compared  to  the  real  part,  k(Tj^)  =  co/c(Tj^),  and  (3p/3x)/p  =  ± 
ico/c(Tj^).  Under  these  conditions  Eq.  (10)  becomes 

ypm\  ^  ^  3x  [  '  - 1)  fx/  (1  +  i)a/5K .  J  (1  +  i)a/5K  / 

Finally,  referring  to  Appendix  B,  we  can  compute  k  and  z  according  to 

=  co^  C(co)  Pc  ,  (12) 

and 

z  =  (k/o))/f2C(0).  (13) 

Note  that  now  C(co),  and  therefore  k,  depends  upon  the  direction  of  propagation  and  that  where  the 
±  appears,  so  does  dTj^j/dx.  The  term  dT^/dx  is  positive  when  the  temperature  increases  in  the 
direction  of  propagation. 

The  presence  of  a  temperature  gradient  establishes  a  difference  in  propagation  constants 
beyond  a  simple  change  in  direction.  In  Region  3, 

pj(3)  =  p.(3)e  i(ot-k3x)  (14) 
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In  the  direction  of  propagation  of  the  incident  wave,  dTjjj/dx  is  positive.  The  reflected  wave  is 

p/3)  =  p^(3)e  i(at-k’3x).  (j5) 

For  dT^/dx  =  0,  k3=k'3  =  0)/c  -  ia.  In  the  presence  of  a  temperature  gradient,  k3  will  have  an 
increasing  imaginary  part  while  k'3  will  have  an  imaginary  part  that  decreases  to  zero  then  becomes 
positive. 

The  only  remaining  task  is  to  get  an  expression  for  impedance  z  is  only 

associated  with  a  boundary  condition,  then  dTj^/dx  should  not  enter  the  calculation  of  z,  i.e.,  z  is 
associated  with  a  boundary,  not  with  a  wave.  This  means  that  at  x  =  I3,  Z3  =z'3  =  z  calculated  for 
that  temperature.  At  x=l4,  Z4  =  z’4  =  z  calculated  at  Tj^=Tho^.  To  include  this  asymmetry,  we 
must  rewrite  the  boundary  conditions  at  x  =  I3  as 

p.(3)e  -ik3l3  +  p^(3)  e  ik'3l3  =  p.(4)  g  -il:4l3  +  e  (16) 

and 

X  p^)  e-ik3l3  .  X  pp)  eik3l3  =  X  [pW  g-Xb  -  eik4l3],  (17) 

or  rewriting. 


and 


pf  ^  e-iJ^sb  -  B  g-  e-iksb]  =  -  P^^  eij^sb  +  B  ^  eiksb 


p(3)_e-ik3l3  [B-1] 

^  ~  eik^b 


(18) 

(19) 


where  B  is  defined  in  Appendix  A.  Working  toward  the  driving  end,  we  must  apply  a  similar 
treatment  at  x  =  I2, 

pW  g-ikjj  +  pp)  gikj2  =  pp>  e-W2  +  pp)  gik;i2  (20) 

and 


Xjp(2)g.ikJj  p(2)gik3J2)  =  Xpp)g-ik5l2.X  pp)giJc;i2  (21) 

‘  ^  Z2  z'3 
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or 


ei^2l2  =  Z2  e-i^sh  -  4-  e-i^sh-i^sb  ei^3^2 
‘  z'-5  z,  LB+IJ 


(23) 


Div^■ding,  and  setting  the  ratio  on  the  right  hand  side  equal  to  C, 


or, 


or 


Pp>  e-i*^2l2  +  Pp)  ei^^2l2  =  CP^  e-i^2l2  -  CPp^  eiJ^2l2 


(24) 


pp)  e-i^2l2  [1-C]  =  -PP^  t^ih  [1+C] 


pp)  =  r^U.]  e-2ik2l2  Pp>. 
LC  +  IJ 


(25) 

(26) 


The  partem  for  no  temperature  gradient  derived  in  Appendi)^  A  is  now  reestablished  and  P^qj^i  has 
the  same  form  as  in  Appendix  A. 

For  the  purpose  of  the  calculations,  the  speed  of  sound,  viscosity,  etc.  within  the  stack  are 
computed  at  the  average  temperature  of  the  stack.  The  results  of  the  calculations  are  shown  in  Figs 
4-8,  along  with  the  experimental  data.  These  figures  are  graphs  of  1/Q,  which  is  proportional  to 
the  net  attenuation  in  the  resonator,  versus  temperature  difference,  ranging  from  zero  to  that 
required  for  onset  of  self  oscillation.  The  figures  correspond  to  mean  gas  pressures  of  170,  238, 
307,  376,  and  500  kPa,  respectively.  The  individual  points  represent  the  measurements  and  the 
lines  show  the  results  of  the  calculations.  Although  the  agreement  between  the  theory  and 
experiment  is  generally  pleasing,  there  are  some  noticeable  differences.  The  theory  over  predicts 
the  attenuation  at  low  temperature  differences,  though  this  discrepancy  diminishes  with  increasing 
mean  pressure.  Also,  the  predicted  onset  temperature  difference  increases  with  mean  pressure, 
ranging  from  approximately  320  K  at  170  kPa  mean  pressure  to  460  K  at  500  kPa.  The  measured 
onset  temperature  difference  does  not  show  such  a  regular  pattern. 

Measurements  were  also  made  for  the  second  and  third  longitudinal  modes  of  the  resonator  at 
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The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  170  kPa. 
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t,urc  6.  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover  stack.  The 
symbols  represent  tlie  data,  while  the  line  represents  the  results  of  the  calculations. 
The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  307  kPa. 


O/l- 


67 


0.02 


0.01 


ic  prune  mover  is  filled  with  helium  at  a  mean  pressure  of  376  kPa. 


0.05 


O/l 


68 


Figure  8;  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover  stack.  The 
symbols  represent  the  data,  while  the  line  represents  the  results  of  the  calculations. 
The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  500  kPa. 


two  mean  gas  pressures,  170  and  307  kPa.  The  data  and  calculations  arc  shown  in  Figs.  9  and  10. 
The  X's  and  solid  lines  correspond  to  the  second  mode,  while  the  open  squares  and  dashed  lines 
correspond  to  the  third.  It  was  not  possible  to  make  measurements  near  onset  for  these  modes, 
because  the  onset  temperatures  exceed  that  of  the  fundamental.  The  agreement  is  good.  The 
tendency  to  over  predict  the  attenuation  at  zero  gradient  decreases,  especially  for  the  third  mode. 
No  conclusions  can  be  drawn  concerning  the  ability  to  predict  the  onset  temperature. 

III.  SUMMARY  AND  DISCUSSION 

The  frequency  response  of  a  thermoacoustic  prime  mover  has  been  measured  as  a  function  of 
the  mean  gas  pressure  and  temperature  gradient  across  the  prime  mover  stack.  The  frequency 
response  of  two  simpler  resonator  configurations  was  also  measured  in  the  absence  of  an  applied 
gradient.  The  quality  factor  Q  and  resonance  frequency  were  determined  from  the  response. 
Rather  than  using  a  standing  wave  analysis  such  as  that  given  by  Swift,^  we  have  analyzed  the 
results  in  terms  of  counterpropagadng  plane  waves,  an  approach  used  in  studies  of  the  acoustic 
properties  of  porous  materials.  The  motivation  behind  this  analysis  is  to  develop  the  ability  to 
transfer  extensive  research  of  porous  materials  to  performance  predictions  of  thermoacoustic 
devices.  This  type  of  analysis  will  be  especially  useful  when  considering  the  use  of  alternative 
geometries  in  the  prime  mover  stack. 

In  general,  the  predictions  of  the  counterpropagating  plane  wave  analysis  are  in  good 
agreement  with  experiment.  The  predictions  of  Q  and  the  change  in  resonance  frequency  with 
mean  gas  pressure  arc  within  approximately  5%  and  0.4%  of  measured  values  for  the  no 
temperature  gradient  cases.  In  the  cases  where  temperature  gradients  arc  present,  the  agreement  is 
pleasing,  but  there  are  some  noticeable  discrepancies.  Surprisingly,  some  of  the  worst  agreement 
is  for  small  temperature  gradients.  The  reason  for  this  discrepancy  is  unknown.  However,  the 
effect  is  less  noticeable  at  higher  mean  gas  pressures  and  frequencies.  The  analysis  provtides  the 
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Figure  9:  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover  stack  for  the 
second  and  third  longitudinal  modes.  The  X's  and  solid  line  correspond  to  the 
second  mode,  while  the  open  squares  and  dashed  line  correspond  to  the  tliird.  The 
prime  mover  is  filled  with  helium  at  a  mean  pressure  of  170  kPa. 
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Figure  10;  Graph  of  I/Q  vs  the  temperature  difference  across  the  prime  mover  stack  for  the 
second  and  third  longitudinal  modes.  The  X's  and  solid  line  correspond  to  the 
second  mode,  while  tlie  open  squares  and  dashed  line  correspond  to  tlie  tiiird.  The 
prime  mover  is  filled  with  helium  at  a  mean  pressure  of  307  kPa. 


best  overall  agreement  at  intermediate  gas  pressures.  There  is  a  tendency  to  under  predict  the  onset 
temperature  at  lower  pressures  and  over  predict  it  at  higher  pressures.  We  have  not  presented  data 
or  predicdons  for  the  resonance  frequency  as  a  function  of  temperature  gradient.  However,  the 
predictions  agree  with  measured  values  to  within  a  few  Hz,  a  discrepancy  of  approximately  1  %. 

Even  taking  into  consideradon  these  discrepancies,  our  analysis  provides  a  good  overall 
predicdon  of  the  performance  of  the  prime  mover  below  onset.  However,  the  true  test  is  predicdng 
the  performance  above  onset.  These  are  several  significant  differences  between  prime  movers 
above  and  below  onset.  Most  obvious  is  that  the  acousdc  amplitudes  above  onset  are  quite  large. 
Rados  of  acoustic  pressure  amplitude  to  mean  gas  pressure  of  1  -  10  %  are  common.  Yet,  our 
analysis  is  based  on  linear  acousdcs.  Beyond  the  theoretical  complications  added  by  introducing 
nonlinear  effects,  there  are  a  number  of  experimental  complications.  For  instance,  the  acoustic 
displacement  amplitudes  approach  the  length  of  the  heat  exchangers.  Also  because  of  the  large 
transport  of  heat,  there  is  no  assurance  that  the  temperature  distribution  along  the  prime  mover 
stack  will  be  uniform. 

The  next  logical  step  in  this  research  is  to  investigate  prime  movers  above  onset,  but  still  in 
the  linear  acoustics  regime.  This  study  requires  the  use  of  some  type  of  known  additional 
attenuation  to  limit  the  amplitudes.  Also,  the  counterpropagating  plane  wave  analysis  should  be 
compared  to  Swift's  analysis.^ 
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APPENDIX  A:  Derivation  of  P, 


total 

The  derivation  of  Pj^^g^  in  the  absence  of  a  temeprature  gradient  is  presented  in  this  appendix. 
The  combination  of  elements  in  the  thermoacoustic  prime  mover  creates  a  series  of  boundaries, 
each  with  a  characteristic  impedance.  The  total  acoustic  field  inside  the  resonant  tube  will  be 
evaluated  by  considering  the  field  as  a  superposition  of  traveling  waves  moving  in  each  direction 
through  this  series  of  boundaries  and  porous  elements. 

The  incident  and  reflected  waves  in  the  five  regions  can  be  expressed  as 

pp)  =  pi(i)  e>X“t-kix)  p^(i)  =  p^(i)  girc.+k-ix) 


p-(5)  =  p-(5)  eXmt-kjx)  p^(5)  =  p^(5)  ei(co:+k5x) 

First  consider  the  rigid  end  at  x  =  I5.  For  a  rigid  boundary,  u=0.  Newton’s  law  gives 
po  du/dt  -  -V  p.  Assuming  time  harmonic  plane  waves,  3u/3t  =  jcou  so  jcopgu  =  -V  p.  Hence, 
the  boundary  condition  u  =  0  requires  that  V  p  =  0  at  x  =  I5  or 

[V  (PX^)  +  p/5)  eX©t+k5x))  ]  =  0  ,  (Al) 

X=l3 

or 

P/5)  =  e-2*sh  (A2) 


Next  consider  the  boundary  at  x  =  I4.  Continuity  of  pressure  gives 

p.(4)  ei(ot-k4l.)  +  p^(4)  £i(cat+kJU)  =  pj(5)  ^iCw-kjL.)  +  Pj,(5)  eiCut+kjh)^ 

or  using  Eq.  (A2) 


(A3) 
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pX4)  e-iluU  +  p^(4)  gikJU  z=  Pj(5)  [g-iksU  +  g-SiksIs  + 


(A4) 


Continuity  of  volume  velocity  needs  to  take  into  account  changes  in  porosity  and  can  be  written  as 

[Ui(4)+Ur(4)  =  Ui(5)  +  Ur(5)]  _  (A5) 

But,  U  =  Su  =  Sp/Zgg^,  where  S  is  the  cross  sectional  area  of  the  fluid  and  is  the  specific 
acoustic  impedance  of  the  region.  Hence,  continuity  of  volume  velocity  can  be  written  as 


S^p.(4)  .  P^(4)  g+iktl,)  =  S^p^(5)  g-iJcjU  .  p^(5)  e+ikjU) 


7(4) 


t(5) 

^sa 


(A6) 


Dividing  by  the  cross  sectional  area  of  the  empty  tube  S,  defining  the  porosity  =  S^”VS,  and 
further  defining  the  impedance  z  =  yields 


X  [pj(4)  e-iiuU  .  p/4)  P^[e-ik5U  .  e-2iks!5 


(A7) 


Dividing  Eq.  (A4)  by  Eq.  (A7)  and  solving  for  gives 


where 


p/4)  =  e-2ik«U  J  p.(4)^ 


^  _  Z5  +  e~2W5 

Z4  e-ikjl»  _  g-ZiWj  eikjU 


(A8) 


(A9) 


Now  to  the  boundary  at  x  =  I3; 

p.(3)  e-ik3l3  +  Pj(3)  eik3l3  =  Pj(4)  e-ikJ,  +  p^(4)  ei}-j3,  (A  10) 

and 

X[pj(3)  e-ik3]3  .  p/3)  5,i}:3lj]  ^  X[p.(4)  e-ibls  .  p/4)  giJalj],  (AH) 

Z3  Z4 

so 
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(All) 


1 

+  Q-2iktU 

[A-r 

-B-l 

A+l. 

Z3 

j  Q-iktl}  _ 

1 

■A-1 

■A+lJ 

e^h| 

=  B 


and 


p(3)  =  p.(3), 

"  Ib+iJ  ‘ 


We  can  now  identify  a  pattern  and  write 


where 


and 


where 


LC+lJ 


C=  ^ 


Z2 


g-ikjli  ^  £-2ik3l3 

[3-r 

Ib+I- 

g-iksU  _  g-2ikslj 

■JzLl 

-B+lJ 

eiksh  1 

p,(i)  =  p,(l), 


D  = 


Z2 


(  g-iJcilt  +  Q-2ikj2 

\£zi: 

Ic+l. 

giJcJi 

,  ^-2i3c2l2 

'£zL] 

■C+lJ 

j 

(A13) 

(A14) 

(A15) 

(A16) 

(A  17) 


Now  we  must  consider  the  driver.  We  will  assume  that  the  driver  delivers  a  constant 
displacement  /  independent  of  load,  x=/  e*^^^  In  this  case,  the  velocity  (provided  I  is  small)  is 
Uo(Oit)  =  "^(^1  Applying  continuity  of  volume  velocity  (with  the  definition  of  Zj^)  at  x=0. 


so 


and 


Ptoui  =  +  P/^^  =  iO)/zi 


1  +  e’^^>^' 

[D-l] 

lD+1 

1  - 

D-l] 

•D+lJ 

(A18) 

(A  19) 


(A20) 
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APPENDIX  B:  Parallel  Slits  Without  a  Temperature  Gradient 


Using  fundamental  concepts  from  inviscid  acoustics,  it  is  easy  to  show  that  the  propagation 
constant^and  specific  acoustic  impedance^an  be  expressed  as 


and 


=  0)2  C  pf. 


Zsa  =  (k/0))/C. 


(Bl) 


(B2) 


In  these  equations  o)  is  the  angular  frequency,  pf  the  mean  fluid  density  and  C  the  compressibility 

r 

of  the  medium.  One  way  of  defining  the  compressibility  is  through  the  equation  of  state  p  =  s/C,  p 
is  the  acoustic)^j5ressure  and  s  the  condensation.  These  results  must,  of  course,  be  modified  when 
considering  propagation  through  porous  media  to  take  into  acoustic  the  complications  introduced 
by  viscosity  and  thermal  conductivity.  The  approach  commonly  taken  in  porous  media  analysis  is 
to  incorporate  the  added  complexity  into  the  definitions  of  new  parameters,  such  as  the  complex 
density  and  compressibility,  leaving  the  forms  of  Eqs.  (Bl)  and  (B2)  unchanged.  This  analysis  is 
outlined  in  this  Appendix.  The  reader  is  directed  to  Refs.  7,  9  and  10  for  a  more  complete 
analysis. 

The  average  acoustic  velocity  in  the  slit  formed  between  two  stationary  parallel  plates 
separated  by  2a  is^*^ 

i)a/5v]l 


icopf^x 


(1  +  i)a/5v  } 


(B3) 


where  where  5v  =  (2v/co)^^,  and  v  is  the  kinematic  viscosity.  We  define  the  complex  density  pj. 


in  terms  of  the  pressure  gradient  and  the  slit-averaged  acoustic  velocity  such  that 


^  =  Pcl^|. 


at 


(B4) 


To  be  consistent  we  should  use  the  slit- averaged  pressure  gradient  in  Eq.  (B4).  However,  we 
have  assumed  the  pressure  to  dependent  only  on  x,  so  this  distinction  need  not  be  made.  Solving 
Eq.  (B3)  for  -{dp/dx)  and  comparing  the  result  to  Eq.  (B4)  gives 
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/,  taiih[(l  +  i)a/5v 

Pc  =  Pf  1 - — 

\  (1  +  i)a/5v 

The  temperature  equation  is 

3^0 

icopfCpS  +  icop  =  Xh - (B6) 

where  is  the  coefficient  of  thermal  conduction,  0  the  difference  between  slit  wall  and  fluid 
temperatures,  and  Cp  the  specific  heat  of  the  fluid,  y  is  measured  perpendicular  to  the  slit  wall. 
Applying  the  boundary  condition  that  0  =  0  at  a  yields  the  solution 


0  =  p/pfCp 


^  cosh[(l +i)y/5j^ 
cos  h[(l+i)a/5jj 


(B7) 


A  differential  form  of  the  ideal  gas  equation  of  state  is  dp/pf  =  dp/p^^^  -  dT/Tjj,_  Substituting 
acoustic  quantities  for  differentials  and  identifying  the  condensation  as  p/pp  we  have 

s  =  p/p^-drr^.  .  (B8) 


Substituting  Eq.  (B7)  into  Eq.  (B8),  and  noting  that  Pj^p^CpT^  =  (y-  l)/Yfor  an  ideal  gas,  gives 


s  = 


P 

TPm 


[ 

l  +  (y-l) 


cosh[(l  +  i)y/5K]  j 
cosh[(l  +  i)a/5j/ 


Averaging  the  condensation  across  the  slit  gives. 


TPm  \  [(1  +  i)a/5j  / 


(B9) 


(BIO) 
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Zwikker  and  Kosten^  define  the  complex  compressibility  of  the  gas  in  the  slit  through  the 
silt-averaged  equation  of  state 

P  =  s/C(£0)  (Bll) 

Substituting  Eq.  (BIO''  into  Eq.  (Bll)  and  keeping  in  mind  that  the  pressure  is  assumed  constant 
across  the  slit,  we  arrive  at  the  following  expression  for  the  complex  compressibility 

=  (B12) 

[a+i)a/5j  I 

Now  that  we  have  the  (slit-averaged)  complex  density  and  compressibility,  we  can 
immediately  write  down  the  slit-averaged  expressions  for  the  propagation  constant  and  impedance. 
From  Eq.  (31), 


k^  =  co2  _  (3J3) 

From  Eq.  (B2)  and  the  relation  z  =  (see  Appendix  A) 

z  =  (k/co)/nC(a)).  (B14) 

To  find  common  ground  with  those  readers  familiar  with  Swift’s  work,  especially  Ref.  1, 
our  Eqs.  (B3),  (B6),  and  (B7)  are  identical  to  (from  Ref.  1)  the  slit  averaged  Eq.  (A4),  and  Eqs. 
(A9)  and  (A  10)  in  the  absence  of  a  temperature  gradient,  respectively. 
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Argon 


Pressure  (kPa) 

f^Cmeasured) 

foCprcdicted) 

%  error 

Q'meas'jred) 

Q(predjci£d) 

%  error 

170 

155.6 

156. 

+0.3 

79 

81 

+2.5 

238 

155.8 

156.5 

+0.4 

91 

95 

+4.4 

307 

156 

156.5 

+0.3 

102 

no 

+7.8 

376 

156,1 

+QJ 

136 

123 

+6.0 

Af  =  0.5  Af  =  0.5  avg  =  +0.3  avg  =  +5.l 


Helium 


Pressure  (kPa) 

fQ(measured) 

fQtjredicled) 

50  error 

Q(measuied) 

CXpredicied) 

%  error 

170 

489 

437.5 

+0.3 

48 

50 

+4.2 

238 

490.2 

488 

+0.4 

57 

59 

+3.5 

307 

490.6 

488.5 

+0.4 

64 

67 

+4.7 

376 

490,5 

482 

iSLl 

70 

74 

±12 

Af=  1.6 

Af=  1.5 

avg  =  +0.35 

avg 

=  +4.5 

Table  I.  Measured  and  predicted  values  of  the  resonance  frequency  and  Q  the  empty 
resonator  for  helium  and  argon  gas  at  various  pressures.  Af  is  the  total  change  in 
resonance  frequency  with  pressure. 
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Argon 


Pressure  (kPa) 

fo(measured) 

fo(predicted) 

%  error 

Q(measuied) 

Qipredicted) 

%  error 

170 

156.4 

156.5 

+0.1 

60 

59 

-1.7 

238 

156.5 

157 

+0.3 

72 

71 

-1.4 

307 

156.4 

157 

+0.4 

81 

82 

+  1.2 

376 

156.7 

157 

+0.2 

89 

91 

+2.2 

500 

157.1 

m 

iQJ. 

102 

103 

±Ln 

11 

P 

o 

II 

< 

avg  =  +0.2 

avg 

=  +0.3 

Helium 


Pressure  (kPa) 

f^Cmeasured) 

fQ(predicted) 

terror 

Q(measured) 

CXpredicted) 

%  error 

170 

490.3 

488.5 

-0.4 

36 

35 

-2.8 

238 

491.4 

489.5 

-0.4 

43 

42 

-2.3 

307 

492.0 

490 

-0.4 

49 

49 

0 

376 

mi 

Af=2.8 

mj. 

Af=2.0 

HI 

avg  =  -0.4 

55 

54 

< 

VQ 

11 

ufe 

Table  II.  Measured  and  predicted  values  of  the  resonance  frequency  and  Q  of  the  resonator 
with  the  ambient  heat  exchanger  for  helium  and  argon  gas  at  various  pressures.  Af  is  the 
total  change  in  resonance  frequency  with  pressure. 
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Argon 


Pressure  (k-Pa) 

fo(rrieasured) 

/^(predicted) 

%  error 

QCmeasured) 

Q(predicted) 

%  error 

170 

161.6 

162.0 

+0.2 

40 

38 

-5 

238 

161.8 

162.0 

+0.1 

47 

47 

0 

307 

162.1 

162.0 

+0.0 

54 

53 

-2 

376 

162.1 

162.5 

+0.2 

58 

59 

+2 

500 

162.3 

mi 

±(li 

68 

70 

±2 

Af=0.7  Af  =  0.5  avg=:+0.15  avg  = +2 


Helium 

Pressure  (kPa) 

f^Cmeasured) 

fQ(predlcted) 

%  error 

Q(measured) 

Q{predjcted) 

%  error 

170 

505.2 

504.5 

-0.1 

23 

20 

-13 

238 

507.2 

505.5 

-0.3 

28 

26 

-7 

307 

508.5 

506 

-0.5 

32 

30 

-6 

376 

509.4 

507 

-0.5 

36 

34 

+5 

500 

mi 

Af=5.3 

mi 

Af=3.0 

avg  =  -0.4 

42 

40 

avg  =  -7 

Table  rH.  Measured  and  predicted  values  of  the  resonance  frequency  and  Q  of  the  prime 
mover  with  no  temperature  gradient  for  helium  and  argon  gas  at  various  pressures.  Af  is 
the  total  change  in  resonance  frequency  with  pressure. 
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tanh  (l-i-i)  a/Sy] 

(1+i)  a/5y 

tanh[(l+i)a/SK] 
tanh[(l+i)a/5s] 

5v  =  V2v7  CO 

5s  =  y2Ks/CL) 

=  ^h/PfCp 

Xj  =  solid  thermal  conductivity 
Cj  =  solid  specific  heat/unit  mass 


tanh  [  (1+i)  a/S)j 
(l+i)a/5j: 


5jt  =  ilKf/ca 
~  ^s/Ps'-s 

=  thermal  conductivity  of  the  fluid 

Ps  =  solid  density 
b  =  half  width  of  solid 


Table  IV.  Symbols  used  in  Eq.  (6). 
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ABSTRACT 

An  experimental  and  theoretical  investigation  of  sound  propagation  in  a  porous 
sample  composed  of  capillary  tubes  with  rectangular  cross  sections  is  described  in  this 
paper.  An  experimental  technique  valid  for  low  flow  resistivity  and  high  porosity  porous 
samples  was  developed  to  measure  the  attenuation  and  phase  velocity  in  the  porous 
material.  This  technique  uses  transmission  of  a  short  pulse  in  a  large  tube  through  the 
porous  sample  and  subsequent  frequency  domain  analysis  in  the  range  200  to  1300  Hz. 
Good  agreement  was  obtained  if  an  anomalous  tortuosity  factor  of  1.1  is  used  in  the 
theory.  A  scaling  factor  for  relating  cylindrical  and  square  tube  capillary  theories,  known 
as  the  dynamic  shape  factor,  was  investigated.  Propagation  constants  computed  from  use 
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of  a  near  unity  dynamic  shape  factor  in  the  cylindrical  pore  theory  agree  favorably  with 
calculations  based  on  the  square  pore  theory  for  the  frequencies  and  pore  radii  used  in  the 
experiment. 

(PACS  #’s:  43.28.Fp,  43,20.Bi,  43.20.Mv,  43.85.Dj) 

INTRODUCTION 

The  interaction  of  sound  with  porous  media  has  many  practical  applications  and  has  a 
long  history.  Porous  media  are  ideal  sound  absorbers  for  use  in  architectural  acoustics  and 
anechoic  chambers.  The  porous  nature  of  the  Earth's  surface  and  ocean  bottoms  can 
greatly  influence  the  propagation  of  sound  in  the  air  and  oceans.  An  example  of  the  diverse 
uses  of  sound  propagation  in  porous  media  comes  from  our  laboratory  where  we  recently 
used  measurements  and  theory  for  sound  propagation  in  the  porous  ground  to  determine 
physical  parameters  which  are  related  to  the  agricultural  suitability  of  the  soils 
investigated.^ 

Several  models^'^  for  porous  media  are  based  on  the  adaptation  of  the  solution  for 
sound  propagation  in  cylindrical  capillary  tubes  to  pores  of  irregular  geometries.  The  basic 
fluid  field  equations  used  in  these  models  are  a  simplified  version  of  the  linearized  Navier- 
Stokes  equations^  for  a  fluid.  The  approximation  employed  in  this  calculation  is  that  the 
transverse  fluid  velocities  are  much  smaller  than  the  longitudinal  fluid  velocity.  Zwikker 
and  Kosten^  were  first  to  obtain  the  solution  based  on  the  simplified  version  of  the  fluid 
model  equations  and  showed  that  these  solutions  agreed  with  Kirchoffs  exact  solution‘s  in 
the  limit  of  high  and  low  frequencies.  Tijdeman^  investigated  the  range  of  validity  of  the 
Zwikker  and  Kosten  solution  in  comparison  to  the  more  rigorous  Kirchoff  solution  and 
determined  that  the  condition  for  the  approximation  to  hold  was  the  condition  on  velocities 
listed  above.  He  refers  to  the  Zwikker  and  Kosten  theory  as  the  "low  reduced  frequency 
approximation”. 
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We  have  developed  a  porous  media  model  for  rectangular  cross  section  capillary 
tubes.  The  rectangular  pore  calculation  is  based  on  the  low  reduced  frequency 
approximation.  Specifically,  the  model  allows  one  to  compute  propagation  constants  and 
characteristic  impedances  for  porous  media  consisting  of  rectangular  pores.  The  effects  of 
a  slight  capillary  tube  curvature  or  tortuosity  are  accounted  for  in  this  model.  For  the  single 
pore,  a  series  solution  is  used  for  the  particle  velocity,  pressure,  density,  and  temperature. 
The  extension  of  the  single  pore  theory  to  bulk  media  consisting  of  an  array  of  pores  is 
developed  using  standard  techniques. This  is  discussed  in  Sec.  I.  The  single  pore 
solution  for  rectangular  pores  has  been  developed  independently  by  Stinson.^ 

The  rectangular  pore  model  is  compared  to  attenuation  and  phase  velocity 
measurements.  The  ceramic  porous  samples  used  in  the  experiment  consisted  of  nominally 
straight  capillary  tubes  having  square  cross  sections.  These  low  flow  resistivity,  high 
porosity  samples  may  be  useful  for  low  frequency  sound  absorption.  A  more  thorough 
description  of  the  po.  ous  samples  and  experimental  technique  is  given  in  Sec.  11. 
Comparison  of  theory  and  experiment  is  discussed  in  Sec.  III. 

A  general  theory  for  arbitrary  pore  shape^  was  developed  by  introducing  a  scale 
factor  known  as  the  dynamic  shape  factor  to  scale  between  different  pore  geometries.  The 
limiting  cases  for  this  scaling  factor  are  circular  pores  and  parallel  slits.  The  limiting  cases 
of  rectangular  pores  are  square  pores  and  parallel  slits,  which  give  rectangular  pore  theory  a 
wide  range  of  applicability.  The  dynamic  shape  factor  for  square  pores  is  frequency 
dependent  just  as  it  is  for  parallel  slits. ^  A  discussion  of  the  dynamic  shape  factor  for 
square  pores  is  given  in  Sec.  IV. 

I.  PROPAGATION  IN  RECTANGULAR  PORE  .MEDIA 

Sound  propagation  in  a  single  rectangular  capillary  tube  is  developed  first.  The 
acoustic  field  in  the  pore  is  specified  to  first  order  in  the  acoustic  variables.  Boundary 
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conditions  at  the  pore  wall  are  that  the  walls  are  rigid  and  thus  the  total  particle  velocity  is 
taken  as  zero.  Due  to  the  high  heat  capacity  and  thermal  conductivity  of  the  pore  wall,  the 
temperature  of  the  fluid  in  the  pore  at  the  boundary  is  taken  to  be  the  same  as  the  pore  wall. 
We  do  not  assume  any  internal  mean  flow.^  The  single  pore  theory  is  then  used  to  develop 
the  theory  for  a  porous  material  consisting  of  an  array  of  rectangular  capillary  tubes 
possibly  having  a  slight  longitudinal  curvature  or  tortuosity. 

A.  Acoustical  disturbances  in  a  single  rectangular  capillary  tube. 

The  coordinate  system  shown  in  Fig.  1  has  the  z  axis  parallel  to  the  tube  axis.  The 
transverse  dimension  is  spanned  by  an  x-y  coordinate  system  with  origin  at  the  lower  left 
comer  of  the  rectangle.  First  order  acoustic  variables  are  the  real  parts  of: 


p(z,t)  =  Po  +  Pi(z)  exp(-icot). 

(1) 

v(x,y,z,t)  =  [Vj^(x,y,z),  Vy(x,y,z),  v^(x,y,z)]  exp(-icot) 

(2) 

T(x,y,z,t)  =  Tq  +  Tj(x,y,z)  exp(-ia)t). 

(3) 

and 


p(x,y,z,t)  =  Po  +  Pi(x,y,z)  exp(-iwt).  (4) 

Subscript  0  refers  to  ambient  values.  Subscript  1  implies  first  order  and  in  Eq.  (2)  the  x,y, 
and  z  components  of  particle  velocity  are  v^,  v^,  and  v^,  respectively.  Equations  (1),  (3), 
and  (4)  are  the  acoustic  pressure,  temperature,  and  density.  Acoustical  disturbances  within 
the  pore  are  taken  to  satisfy  the  following  relations: 
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-iwPoV^(x,y,z)  =  -  +  Ti 


\dx^  dy^j 


v,(x,y,z). 


(5) 


-iwpi(x,y,z)  +  pQ 


^3vx(x,y,z)  ^  3vy(x,y,z)  ^  8vz(x,y,z)^ 


V 


8x  3y 


3z 


0 


(6) 


y 


Y 


pj(x,y,z)  =  -  pQpTi(x,y,z)  +  — Pj(z), 


(7) 


and 


-io^p^c  Tj(x,y,z)  =  -icopT^PjCz) 


+  K 


8x2  3y2 


Tj(x,y,z). 


(8) 


These  relations  are  given  in  the  frequency  domain  where  8/8t  is  replaced  by  -ico  where  t  is 
time  and  O)  is  radian  frequency.  Response  functions  and  transport  coefficients  are  Cp,  the 
constant  pressure  heat  capacity  per  unit  mass;  y,  the  ratio  of  specific  heats;  c,  the  adiabatic 
sound  speed;  p  =  -  (8p/8T)p/pQ,  the  thermal  expansion  coefficient;  r\,  the  viscosity;  and  K, 
the  thermal  conductivity.  In  order,  these  equations  express  the  z  component  of  the  equation 
of  motion,  continuity  or  mass  conservation,  equation  of  state  for  density,  and  heat  transfer. 
Equations  (5)  through  (8)  are  the  same  set  of  equations  used  by  Zwikker  and  Kosten.^  In 
using  these  approximate  equations  the  assumption  is  that  the  transverse  velocity 
components  v^  and  Vy  are  much  less  than  the  longitudinal  velocity  v^.  Further  discussions 
of  these  approximations  can  be  found  in  Appendix  B  of  Ref.  6  and  Appendix  A  of  Ref.  10. 

The  following  notation  will  be  used  to  facilitate  comparison  with  Attenborough's 
results. A  dimensionless  "shear-wave  number"  which  is  proportional  to  the  ratio  of  the 
pore  radius  and  the  viscous  boundary  layer  thickness  is  ^  =  Rlp^w/ril^''^.  Here  R  is  a 
characteristic  transverse  dimension  of  the  pore.  For  definiteness  we  take  R  to  be  twice  the 
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transverse  pore  area  divided  by  the  transverse  pore  perimeter.  Thus,  R  is  the  tube  radius 
for  a  cylindrical  pore  and  is  the  semi-width  for  a  square  pore.  "Wide  tubes"  with  the  same 
R  are  acoustically  equivalent,  i.e.  they  have  the  same  propagation  constants.^  A 
dimensionless  number  proportional  to  the  ratio  of  the  pore  radius  to  the  thermal  boundary 
layer  thickness  is  =  RCp^coCp/K)^  ^^  or  where  Npj.  =  riCp/K  is  the  Prandtl 

number. 

To  make  rapid  progress,  denote  the  z-component  of  panicle  velocity,  v^,  by 


Vj,(x,y,z)  = 


F(x,y;:\.)  dpi(z) 
icopo  dz 


Similarly,  denote  by 


(9) 


Ti(x,y,z)  =  -^Pl(z)  F(x,y;XT) 

c2poP 


(10) 


the  acoustic  pore  temperature.  The  thermodynamic  relation  T^P^/Cp  =  (y-  l)/c^  can  be 
used  in  Eq.  (8)  for  Tj.  Particle  velocity  and  temperature  ^  (5)  and  (8)  reduce  to  the 

simple  forms 


F(x,y;?i) 


[B^i  ^  ay2j 


F(x,ya)  =  1  , 


and 


F{x,y,Xj)  + 

iA-j 


3X2  0y2^ 


F(x,y;?i-j-)  =  I  , 
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respectively,  subject  to  the  boundary  condition  F  =  0  at  the  pore  boundary.  The  solution 
for  F(x,y;^)  for  rectangular  pore  boundaries  is^* 


F(x,y;X)=lf  £ 

^  m.nodd 


sin(m7ix/2a)  sin(n7cy/2b) 
mnYm,n(^) 


(13) 


where 


=  1  +  (irrW)  b^m2_+  a^.  (14) 

(a+b)2 

and  the  ratio  of  transverse  pore  area  to  perimeter  is  R  =  2ab/(a  +  b).  (Recall  that  X  = 
R(PqC0/ti)^^^).  For  later  use,  the  average  (l/4ab)  j  F(x,y;^)dxdy  over  the  pore  cross 
section  is 


F(^)  =  ^  X  — ■ — 

^  m.nodd  rn^^m  n(^) 


(15) 


The  z-component  of  particle  velocity  is  given  by  Eq.  (9’)  with  F(x,y;?i)  in  Eq.  (13)  and  the 
excess  temperature  by  Eq.  (10)  with  the  replacement  of  A.  in  Eq.  (13)  by  Xj. 

To  derive  a  wave  equation  for  the  pressure,  the  fluid  equations  in  Eqs.  (5)  through 
(8)  are  averaged  over  the  pore  cross  section.  Denote  by  Pi(z)  =  (l/4ab)  I  Pi(x,y,z)dxdy 
the  average  of  the  acoustic  density  in  the  pore,  and  similar  notation  for  v^(z)  and  TjCz)  for 
the  transverse  area  average  of  vz(x,y,z)  and  T](x,y,z).  Use  of  Eq.  (9)  and  Eq.  (10)  for  the 
z-component  of  particle  velocity  and  temperature,  and  the  fluid  equations  (5)  through  (8) 
results  in  a  set  of  averaged  equations. 


i5:^v,(z)=^  (16) 

F(A.)  ^ 
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-iwpl(z)  +  Po 


dvz(z) 

dz 


=  0, 


(17) 


and 


Pl(z)  = 


(l-Y)F(XT)  +  y 


Pi(z)  • 


(18) 


The  boundary  conditions,  v^(x,y,z)  =  0  and  Vy(x,y,z)  =  0  at  the  boundary,  were  used  in 
obtaining  Eq.  (17)  from  the  continuity  equation,  Eq.  (6).  Also,  Eq.  (10)  for  Tj(z)  was 
used  in  the  equation  of  state  Eq.  (7)  to  obtain  Eq.  (18).  Following  Attenborough^  we 
define  a  complex  density  from  Eq.  ( 1 6)  and  complex  compressibility  from  Eq.  (18). 


P  = 


Po 

pa)  ’ 


(19) 


and 


g  ^  j_  ^  ^  (1  -Y)F(^t)  +  Y 

Po  Pi  poC2 


(20) 


Eliminating  pj  from  continuity  and  state  Eqs.  (17)  and  (18)  and  using  Eqs.  (19)  and  (20) 
we  obtain: 


i03  p  v^(z) 


^Pi(z) 

dz 


(21) 


and 
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dvzCz) 

dz 


-i(oC  pj(z)  =  0 


(22) 


B  .  Extension  to  bulk  media:  Propagation  constants  and  characteristic 
impedance. 

Consider  a  fluid  half  space  overlying  a  porous  half  space  saturated  by  the  same  fluid. 
The  pores  are  taken  to  have  a  rectangular  cross  section  and  let  the  capillary  tube  axis  of  each 
pore  be  at  an  angle  6  with  respect  to  the  surface  normal.  The  tortuosity  q  =  l/cos9  for  such 
a  porous  sample.^’^>^>^^  The  open  volume  divided  by  the  total  volume  is  the  porosity  Q  of 
the  sample.  The  boundary  conditions  are  continuity  of  volume  velocity  (from  mass 
conservation)  and  the  continuity  of  pressure  (from  Newton's  third  law)  at  the  porous 
interface.^’^’-’^  As  fluid  only  flows  into  the  pores,  the  bulk  particle  velocity  in  the 
porous  media  is  v^/  In  order  to  account  for  propagation  in  a  slanted  pore 

(or  other  tortuous  path),  the  differential  dz  in  Eqs.  (21)  and  (22)  is  replaced  with  qdz  where 
q  >  1  is  the  tortuosity Thus,  the  bulk  acoustical  equations  are 


icopq 


Vzb(z) 

Q 


dpi(z)  _  Q 
qdz 


(23) 


and 


dVzb(z) 

Qdz 


ico  C  Pi(z)  =  0. 


(24) 


Differentiating  Eq.  (23)  by  z  and  eliminating  with  Eq.  (24)  give  an  expression 


d^Pl(z) 

dz^ 


+  (0^  p  C  q2pj(z)  =  0 


(25) 


for  the  pressure  in  the  porous  media.  Assuming  pj  exp(ikz)  gives  the  dispersion  relation 
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-  7  „2  J(1-Y)F(>^t)  +  Y 

=  (0^  p  C  =  ^  q2  - - 

c2  F(>.) 


(26) 


for  the  propagation  constant  k.  From  Eq.  (23),  the  characteristic  impedance  is 


Z  = 


pcoq  _  Po  q 


Qk  F(X)^'^^^  (1-Y)F(Xt)  +  Y 


1/2 


(27) 


The  plane  wave  pressure  reflection  coefficient  Tp  for  a  wave  normally  incident  from  the 
fluid  half  space  on  the  porous  sample  is 


z  -  poc  i 

[m) 

a-Y)F(?tT)  +  Y 

p  Z  +  poC  J 

[f(X)[(1-7)F(Xt)  +  Y 

(28) 


This  relation  will  be  used  in  the  next  section.  Principal  results  of  this  section  are  the 
propagation  constant  Eq.  (26)  and  characteristic  impedance  Eq.  (27)  which  for  rectangular 
pore  porous  media  are  to  be  evaluated  with  the  function  F(X)  given  in  Eq.  (15). 

II.  MEASUREMENT  OF  THE  PHASE  VELOCITY  AND  ATTENUATION 
A.  Description  of  the  porous  media. 

A  schematic  drawing  of  the  porous  sample  is  shown  in  Fig.  2.  Each  subsection  of 
the  composite  sample  was  of  nominal  length  7.68  cm.  Figure  2  indicates  a  composite 
sample  made  from  5  individual  pieces.  Individual  pieces  were  taped  together  at  the  joint. 
A  sheet  of  Teflon  was  wrapped  around  the  composite  to  facilitate  insertion  of  the  sample 
into  the  tube  and  to  seal  the  sample-tube  interface.  Individual  pieces  were  ceramics  made 
by  Coming.  xhe  pores  of  each  piece  were  nominally  square  in  cross  section  and  were 
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nominally  straight  in  the  longitudinal  direction. 

Three  different  square  pore  media  were  investigated.  Table  1  lists  characteristics  of 
each.  In  Table  1,  porosity  was  estimated  using  Q  =  (number  pores/unit  area)  x  (2a)^ 
where  a  is  the  square  semi-width.  Flow  resistivity  was  computed  using  Attenborough's 
equation^  a  =  8qq^S/a^Q  where  S  is  a  steady  flow  shape  factor  S  =  0.89  for  square 
pores and  a  tortuosity  q=l.l  was  estimated  from  fitting  the  theory  and  experiment  for  the 
propagation  constants  as  discussed  below.  In  comparison  with  other  porous  media, ^  the 
square  pore  samples  have  low  flow  resistivity  and  high  porosity.  The  average  semi-widths 
are  listed  in  table  1  for  the  200  and  400  pores/inch^  material  for  which  the  pore  cross- 
sectional  shape  is  well  approximated  as  a  square.  However,  due  apparently  to  differences 
in  the  manufacturing  process,  the  pores  of  the  300  pores/inch^  samples  were  not  well- 
approximated  by  a  square  shape.  Two  opposite  comers  of  the  otherwise  square  shape 
were  rounded.  The  semi-width  of  a  =  0.50  mm  listed  in  Table  1  for  the  300  pores/inch^ 
samples  was  determined  from  the  shortest  diagonal  length  divided  by  ( 2 No  explicit 
use  was  made  of  the  calculated  flow  resistivity  in  the  theory  for  the  propagation  constants. 

B.  Experimental  Apparatus. 

A  block  diagram  for  the  experimental  apparatus  used  to  determine  the  attenuation  and 
phase  velocity  of  sound  in  porous  media  is  shown  in  Fig.  3.  A  single  cycle  of  a  sine  wave 
of  duration  1.3  ms  (for  a  center  frequency  of  750  Hz)  was  generated  using  an  HP  3314A 
Function  Generator.  The  signal  was  amplified  using  a  Krohn-Hite  Model  7500  Amplifier 
and  was  added  to  a  DC  polarizing  voltage  of  125  volts.  The  capacitive  driver  consisted  of 
an  aluminized  mylar  membrane  stretched  over  a  grooved  backplate.^^  A  Teflon  ring 
around  the  perimeter  of  the  driver  was  used  to  hold  the  mylar  in  place  and  to  seal  the  driver 
inside  of  the  tube.  The  tube  was  made  of  aluminum,  had  a  length  of  6.09  meters,  an  inside 
diameter  of  14.6  cm,  and  a  wall  thickness  of  1.11  cm.  Radio  Shack  Model  270-090 
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microphones  were  glued  into  10  cm  long  natural  gas  nipples  having  inside  diameters  of  9 
mm.  Holes  were  made  and  threaded  in  the  tube  60  cm  from  each  end  and  the  microphones 
were  inserted  to  be  flush  with  the  inner  tube  wall.  Microphone  signals  were  amplified 
using  Tektronix  AM  502  Differential  Amplifiers.  A  Masscomp  5550  Minicomputer  with  a 
12  bit  analog-to-digital  board  was  used  to  record  and  analyze  the  microphone  signals.  The 
digitizing  rate  was  300  kHz.  The  HP  Function  Generator  was  used  to  trigger  the 
Masscomp  5550  and  30  pulses  were  averaged  in  the  time  domain  for  each  measurement. 
The  purpose  of  Microphone  2  was  to  give  a  reference  time  and  space  location  for  a  pulse 
traveling  in  the  empty  tube  so  that  the  ambient  sound  speed  could  be  determined  from  the 
pulse  arrival  time  at  Microphone  1. 

A  single  microphone  measurement  method  was  used  to  determine  the  propagation 
constants  in  the  square  pore  media.  Figure  4a  shows  Microphone  1  measurement  of  pulses 
both  with  and  without  a  sample  present  in  the  tube.  As  expected  the  pulse  measured  with 
the  sample  present  is  delayed  in  time  and  attenuated  on  account  of  passage  through  the 
porous  media.  The  Fourier  transform  of  a  typical  pulse  indicated  that  the  pressure  level 
was  about  20  dB  above  the  background  for  frequencies  in  the  range  200  Hz  to  1300  Hz. 
The  cutoff  frequency^  above  which  non-planar  modes  can  propagate  in  the  tube  is  1376 
Hz. 


C.  Determination  of  attenuation  and  phase  velocity  by  the  transfer  function 
method. 

Denote  by  pQ(t)  the  incident  pressure  pulse  at  the  right  end  of  the  sample  and  denote 
by  PpCf)  the  Fourier  transform  of  pQ(t).  The  spectrum  Pj^(f)  after  passage  through  a  sample 
of  length  D^  is 


exp(ikDm) 


1  -  rp2  exp(2ikDm) 


(29) 
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where  is  the  frequency  dependent  pressure  reflection  coefficient  Eq.  (28)  and  k  is  the 
propagation  constant  for  the  porous  media  Eq.  (26).  Since  the  tube  has  a  large  diameter  we 
approximate  the  characteristic  impedance  of  the  tube  by  p^c  where  c  is  the  adiabatic  sound 
speed  in  air.  The  exp(ikDjjj)  factor  accounts  for  propagation  through  the  sample,  the  (1  - 
rp^)  factor  is  for  transmission  into  and  out  of  the  sample,  and  the  denominator  accounts  for 
the  multiple-reflection  of  waves  within  the  sample. 

In  Eq.  (29)  the  subscript  m  refers  to  the  number  of  subsections  used  for  a 
measurement.  For  example,  m  =  5  in  Fig.  2  and  Fig.  3.  We  may  form  a  transfer  function 
hj^(f)  from  the  results  of  two  experiments  on  different  porous  sample  lengths.  The  Fourier 
transform  of  the  time  domain  pulses  recorded  by  microphone  1  in  Fig.  4b  for  two  different 
sample  lengths  gives  a  transfer  function  from  use  of  Eq.  (29), 

hm®  =  ^  -  k„)(D„  -  D,)l  (30) 

where  =  co/c  is  the  wavenumber  for  sound  in  air  and  k  is  given  in  Eq.  (26).  In  Eq.  (30) 
the  single  subsection  spectrum  pi(f)  was  used  as  a  reference  to  divide  out  the  frequency 
response  of  the  capacitive  driver  and  microphone  and  to  divide  out  the  transmission 
coefficient  (1  -  rp^)  in  Eq.  (29).  The  second  part  of  Eq.  (30)  is  an  approximation 
because  we  have  assumed  [(1  -  rp2  explikD^)  /  (1  -  Tp'^explikDj^)]  ~  1.  The  reflection 
coefficient  rp  in  Eq  .  (28)  is  significantly  less  than  1  for  two  reasons.  First  |  F(X)  |  =  1  for 
the  frequency  range  of  the  present  experiment,  as  a  consequence  of  the  low  flow  resistivity 
of  the  square  pore  samples  used.  Second,  for  the  square  pore  material  the  porosity  Q  is 
large.  Hence,  by  Eq.  (28)  we  expect  rp^  =  0  .  Measurements  of  rpwith  Microphone  2  in 
Fig.  3  for  a  nominal  frequency  of  750  Hz  gave  Irp^l  =  (0.03,  0.05,  and  0.03)  for  the  (200, 
300,  and  400)  pores/inch^  samples,  respectively. 
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The  phase  velocity  and  attenuation  are  computed  from  k  =  co/Cpj,  +  ia  and 


Cph(f)  = 


_ 27tf _ 

(Dm  -  Di)'  ^  /m  /n(hm)  +  ko 


(31) 


and 


a(0  =  -  ^  Re 


In  (10) 


(32) 


where  Cp^  is  the  phase  velocity  and  a  is  the  attenuation  in  dB/cm  and  where  Re  /«(h^) 
and  Im  //i(hj^)  refer  to  the  real  and  imaginary  parts  of  /«(hj^),  respectively. 

Equations  (30)-(32)  were  used  to  analyze  the  time  domain  pulses  to  obtain 
experimentally  the  phase  velocity  and  attenuation  for  the  square  pore  media.  In  all  cases, 
we  averaged  over  30  times  domain  pulses  before  taking  transforms.  We  also  used  five 
subsections  so  that  m  =  5  in  Eqs.  (30)-(32)  and  (Dm  -  Di)  =  30.7  cm  nominally.  Since  a 
transfer  function  technique  was  used,  it  was  not  necessary  to  determine  the  pressure 
absolutely  or  to  determine  the  frequency  response  of  the  microphone.  A  central  assumption 
which  was  verified  experimentally  was  the  repeatability  of  any  pulse  measurements  since 
P5(t)  and  pi(t)  were  measured  at  different  times.  The  experimental  results  for  the  three 
different  pore  sizes  given  in  Table  1  are  displayed  in  Fig.  5.  The  experimental  technique 
described  here  is  similar  to  a  method  used  by  Ding'^  to  measure  the  reflection  coefficient  of 
absorbents. 


III.  DISCUSSION  OF  EXPERIMENTAL  AND  THEORETICAL 
ATTENUATION  AND  PHASE  VELOCITY 

The  experimental  and  theoretical  attenuation  and  phase  velocity  were  determined  from 
use  of  Eqs.  (30)-(32)  and  Eq.  (26),  respectively,  and  k  =  (o/Cph  +  ia.  Figures  5a  and  5b 
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are  the  experimental  and  theoretical  attenuation  and  phase  velocity.  For  the  theory,  physical 
constants  used  were  y  =  1.4,  Np  =  0.707,  =  1.2  kg/m^,  t|  =  18.5  pP,  and  the  adiabatic 

sound  speed  c  from  the  propagation  time  of  a  pulse  between  Microphones  1  and  2  in  Fig 
3. 

To  get  the  acceptable  agreement  among  theory  and  experiment  as  indicated  in  Figs.  5a 
and  5b,  a  tortuosity  of  q  =  1.1  was  used.  Referring  to  Eq.  (26)  for  the  propagation 
constant  k,  note  that  the  effect  of  tortuosity,  (which  is  q  >  1),  is  to  increase  the  attenuation 
a  and  decrease  the  phase  velocity  as  one  would  intuitively  expect.  The  calculated 
propagation  constant  for  a  tortuosity  q  =  1  results  in  an  attenuation  10%  lower  than  the 
measured  value  and  a  phase  velocity  10%  above  the  measured  value.  The  discrepancy 
between  experimental  and  theoretical  phase  velocity  in  Fig.  5b  for  the  300  pores/inch2 
samples  having  pores  of  semi-width  0.50  mm  may  be  due  to  the  irregular  shape  of  the  pore 
cross-section  as  described  in  Sec.  HA. 

Pores  which  have  a  tortuosity  other  than  one  have  a  radius  which  is  not  constant 
along  the  pore,  have  a  slight  curvature  or  tilt  with  respect  to  the  axis  normal  to  the  surface, 
or  are  in  a  material  for  which  the  rigid  frame  assumption  is  not  valid. ^  In  the  ceramic 
square  pore  material  used  in  the  experiment,  the  pores  were  straight  and  the  pore  walls  had 
a  density  and  moduli  much  greater  than  that  of  air.  However,  the  semi-width  of  the 
squares  as  determined  with  a  measuring  microscope  at  the  sample  surface  varied  by 
approximately  5%  from  the  average  value.  This  variation  could  result  in  a  tortuosity  other 
than  1  if  the  5%  radius  variation  extended  down  a  single  pore  also.  Another  possibility  for 
a  tortuosity  other  than  one  is  that  the  ceramic  pore  walls  were  also  porous.  Use  of  a 
measuring  magnifying  glass  indicated  that  the  pore  walls  were  indeed  porous,  with  an 
average  pore  diameter  of  about  75  pm  and  pores  as  large  as  100  pm.  Wall  pores  were 
spaced  by  about  150  pm.  Apparently  the  v'all  pores  did  not  connect  adjacent  square  pores. 
Intuitively,  it  seems  that  the  effects  of  porous  walls  would  be  to  increase  the 
compressibility  of  the  gas  and  hence,  decrease  the  phase  velocity.  The  attenuation  w’ould 
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also  increase  which  is  what  was  observed  experimentally.  Champoux^^  reported  a 
tortuosity  of  1.5  for  the  200  and  400  pores/inch^  samples  using  a  non-acoustical 
technique.  It  was  difficult  to  estimate  the  uncertainty  in  this  measurement  since  his 
technique  requires  accurate  knowledge  of  the  sample  porosity  which,  due  to  wall  porosity, 
was  uncertain  for  the  ceramic  samples  used. 

It  was  not  possible  to  align  the  pores  in  neighboring  subsections  of  the  composite 
sample  as  discussed  in  Sec.  IIA.  To  investigate  the  effects  of  misalignment,  attenuation 
and  phase  velocity  were  measured  for  composite  samples  consisting  of  1-5  subsections.  In 
the  single  subsection  measurement  the  reference  for  the  transfer  function  in  Eq.  (30)  was 
taken  to  be  the  empty  tube  signal.  The  discrepancy  between  these  measurements  was  less 
than  3%  and  showed  no  systematic  trends.  Sample  misalignment  was  probably  not  the 
cause  of  an  apparent  tortuosity  greater  than  one. 

IV.  DYNAMIC  SHAPE  FACTOR  FOR  SQUARE  PORES 

The  hypothesis  of  Attenborough's^  cylindrical  capillary-tube  based  porous  m.edia 
theory  was  that  a  circular  pore  of  radius  a/n  could  be  made  acoustically  equivalent  to 
another  pore  of  characteristic  radius  a  by  proper  choice  of  n.  The  condition  for  acoustical 
equivalence  is  taken  to  be 

Im  Pc(X/n)  =It7x  ps(>.)  (33) 

where  recall  that  X  =  aCp^M/ri)^^^  for  a  circular  pore  of  radius  a  or  a  square  pore  of  semi¬ 
width  a  and  subscripts  c  and  s  refer  to  circular  and  square  pores,  respectively.  This 
condition  occurs  since  the  imaginary  part  of  the  complex  density  p  is  very  large  for  small  X 
and  thus  determines  the  behavior  of  the  propagation  constants  and  impedance  for  small  X.^ 
The  complex  density  is  given  generally  for  rectangular  pores  in  Eq.  (19).  For  square 
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pores  of  semi- width  a  the  function  is 


^  m.nodd  m%l' 


1  + 


(m2  +  n2) 


(34) 


from  Eqs.  (14)  and  (15)  for  the  special  case  a  =  b.  For  cylindrical  pores^, 


FJ?i/n)  =1 - ^ 


ji(vrxyn) 


Vr  7Jr\  Jq^VT  ^Vn^ 


(35) 


where  the  dynamic  shape  factor  n  has  been  inserted  and  the  J's  are  Bessel  functions.  The 
range^  of  n  is  thought  to  be  0.5  <  n  <  1 .  The  dependence  of  n  on  the  shear-wave  number  X 
is  shown  in  Fig.  6. 

The  shear-wave  number  range  in  our  experiments  was  4.5  <  ?i  <  18.  As  indicated  in 
Fig.  6,  n  is  a  frequency  dependent  parameter,  which  has  been  previously  noted.^. 
Choosing  a  representative  value  of  >.  =  4.532,  the  corresponding  value  of  the  dynamic 
shape  factor  is  n  =  0.97.  Propagation  constants  are  shown  in  Fig.  7  for  square  and  circular 
pore  theories  where  the  value  n  =  0.97  was  used.  Propagation  constants  for  circular  pores 
were  computed  from  use  of  Eq.  (35)  for  F(Z)  in  Eq.  (26)  for  k.  There  is  less  than  1% 
difference  between  the  propagation  constants  calculated  with  these  theories. 

In  Eq.  (23)  of  Ref.  8  the  author  uses  a  steady  flow  shape  factor  S  to  define  a  flow 
resistivity  for  general  porous  media  as  a  =  8qq^S/a^Q  where  a  is  a  characteristic  dimension 
of  the  pore  equal  to  the  pore  radius  for  circular  pores.  The  DC  mass  flux  through  a  circular 
pore  of  radius  a/S^''^  is  equal  to  that  of  a  pore  of  characteristic  dimension  a,  e.g.  the  semi¬ 
width  a  for  a  square  pore.  If  we  take  the  DC  limit  of  n  in  Fig.  6  (i.e.  X  =  0),  we  get 
Carman’s  value^^  for  the  steady  flow  shape  factor  of  square  pores,  S  =  0.89  =  rij-)(-^''^. 
Hence,  we  are  led  to  conjecture  that  the  steady  flow  shape  factor*  and  the  square  root  of  the 
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dynamic  shape  factor  are  the  same  when  the  dynamic  shape  factor  is  evaluated  at  the  DC 
limit. 

V.  CONCLUSION 

We  have  developed  a  model  for  porous  media  consisting  of  rectangular  pore  capillary 
tubes.  Viscous  and  thermal  dissipation  were  accounted  for.  A  series  solution  was  obtained 
for  the  transverse  variation  of  the  longitudinal  velocity  and  the  excess  temperature.  A 
measurement  technique  using  frequency  domain  analysis  of  short  pulses  propagated 
through  high  porosity,  low  flow  resistivity  samples  was  developed  for  determining 
propagation  constants.  Propagation  constants  were  measured  for  a  ceramic  porous  media 
having  straight  capillary  tube  openings  with  square  cross  sections.  Use  of  an  anomalous 
tortuosity  factor  q  =  1 . 1  resulted  in  favorable  agreement  among  experimental  and  theoretical 
values  of  the  propagation  constants.  It  was  argued  that  the  non-unity  tortuosity  value  was 
due  to  the  finite  porosity  of  the  ceramic  pore  walls.  A  dynamic  shape  factor  n  =  0.97  was 
suggested  as  the  radius  scaling  factor  for  square  pore  and  circular  pore  theories  for  the 
frequency  range  (200  Hz  to  1300  Hz)  and  pore  sizes  (0.50  mm  to  0.77  mm  semi-width 
square  pores)  of  this  investigation.  The  steady  flow  shape  factor  was  found  to  be 
calculable  from  the  DC  limit  of  the  dynamic  shape  factor  for  square  pores.  Future  work 
will  involve  an  investigation  of  the  anomalous  tortuosity  and  propagation  constant 
measurements  on  longer  samples. 
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TABLE  I.  Geometrical  Properties  of  the  Three  Porous  Samples  Used. 


Pores/unit  area 

(inch‘d) 

Semi-Width  a 

(mm) 

Porosity  Q 

(%) 

Flow  Resistivity  o 

(Nm-4  s) 

200 

0.77 

73 

368 

300 

0.50 

47 

1356 

400 

0.57 

81 

606 

105 
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FIGURE  CAPTIONS 

1 .  Coordinate  system  and  geometry  for  the  single  pore  calculation. 

2.  (a)  Side  and  (b)  end  view  of  a  square  po’^e  porous  sample.  To  form  the  total  sample 
several  pieces  of  nominal  length  7.68  cm  were  put  together.  It  was  not  possible  to 
align  the  squares  from  subsection  to  subsection. 

3.  Block  diagram  of  the  apparatus  used  to  determine  the  phase  velocity  and  attenuation 
of  sound  in  low  flow  resistivity,  high  porosity  samples. 

4.  (a)  Sample  waveforms  for  microphone  1  ^>oth  with  and  without  the  porous  material  in 
the  tube.  The  porous  media  was  constructed  as  shown  in  Fig.  2.  (b)  Expanded 
schematic  view  of  the  tube  and  porous  sample  in  Fig.  3.  The  Fourier  transform  after 
the  incident  pulse  passes  through  the  sample  is  Pm(f)- 

5 .  (a)  Experimental  and  theoretical  attenuation  and  (b)  phase  velocity  for  the  3  square 
pore  semi-widths  given  in  Table  I. 

6.  Dynamic  shape  factor  n  to  scale  between  square  and  circular  pores  of  semi-width  a 
and  radius  a/n,  respectively. 

7.  (a)  Attenuation  and  (b)  phase  velocity  computed  using  a  square  pore  of  semi-width  a 
=  0.50  mm  and  a  circular  pore  of  radius  a  =  0.50  mm  /  n  where  n  =  0.97  is  a  dynamic 
shape  factor. 
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SIDE  VIEW: 

5  SUBSECTION- 
COMPOSITE  SAMPLE 


T 

14.2  cm 

i 
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SQUARE  PORE 
COMPOSITE  SAMPLE 


I —  60  cm  — 
1 


CAPACITIVE 

DRIVER 


6.09  meters 


RELATIVE  PRESSURE 
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1000 


1200 


600  800 

FREQUENCY  (Hz) 
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600  800  1000  1200 
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